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O '. Abstract 

We have carried out a new evaluation of the eighth-order contribution to the electron g — 2 
<N ' . 

using FORTRAN codes generated by an automatic code generator gencodeTV. Comparison of the 
"new" result with the "old" one has revealed an inconsistency in the treatment of the infrared 
divergences in the latter. With this error corrected we now have two independent determinations 
of the eighth-order term. This leads to the revised value 1 159 652 182.79 (7.71) x 1CP 12 of the 



electron g— 2, where the uncertainty comes mostly from that of the best non-QED value of the fine 
structure constant a. The new value of a derived from the revised theory and the latest experiment 
is a -1 = 137.035 999 084 (51) [0.37 ppb] , which is about 4.7 ppb smaller than the previous a . 

PACS numbers: 13.40.Em,14.60.Cd,12.20.Ds,06.20.Jr 
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I. INTRODUCTION AND SUMMARY 



The anomalous magnetic moment of the electron has played a central role in testing the 
validity of QED [l, 2|. The test became very stringent when the precision of measurement 
of the electron and positron was improved by three orders of magnitude over the best earlier 
result 



3] by the University of Washington group in the Penning trap experiment [4j 

a e - = 1 159 652 188.4 (4.3) x 1(T 12 [3.7ppb] , 
a e+ = 1 159 652 187.9 (4.3) x 1(T 12 [3.7ppb] , 



(1) 



where a e = (g — 2)/2 and g is the ^-factor of electron. The main source of the remaining 
uncertainty in Eq. ([T]) is the uncontrolled shift of the frequency due to the resonance between 
the electron and the metal cavity of hyperbolic shape. Brown et al. {5] showed that this 
source of uncertainty can be reduced significantly using a metal trap with the cylindrical 
cavity whose resonance structure can be calculated analytically. 

The recent Harvard measurement is based on the cylindrical cavity. Their value an- 



nounced in 2006 is 



6] 



a e (HV06) = 1 159 652 180.85 (0.76) x 10" i2 [0.66ppb] , 



(2) 



which has a 5.5 times smaller uncertainty than the previous measurements listed in Eq. ([[]). 
Very recently, the same Harvard group has succeeded in reducing the uncertainty further by 
a factor 2.7 [Tf] : 



a e (HV08) = 1 159 652 180.73 (0.28) x 10" 12 [0.24ppb] . 



(3) 



To match the precision of the measurement the theory of a e must include radiative cor- 
rections of up to the eighth-order of QED perturbation theory as well as the hadronic and 
weak contributions 



a e = a e (QED) + a e (hadron) + a e (weak). 



(4) 



The hadronic [8j, |9j, lid, [ill] and weak contributions [12| to a e are very small, but not entirely 
negligible relative to the measurement uncertainties (T5]) or (J3j): 



a e (hadron) = 1.682 (20) x 10~ 12 , 
a e (weak) = 0.0297 (5) x 10~ 12 . 



(5) 
(6) 
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The QED contribution a e (QED) can be divided further into four parts taking account of 
the presence of other leptons: 



a e (QED) = A 1 + A 2 {m e /m^) + A 2 {m e /m r ) + A 3 (m e /m^,m e /m T ) 



(7) 



where rrie.m^, and m T are masses of the electron (e), muon (fi) and tau-lepton (r), respec- 
tively. A%, being dimensionless, depends only on the fine structure constant a. A 2 denotes 
contributions from the Feynman diagrams which have closed loops of either muon or tau- 
lepton. A 3 stands for the contributions of the Feynman diagrams which contain both \x loop 
and r loop. Each A t can be calculated by the QED perturbation theory 



A: = A) 



(2) 



a\ 3 

7T- 



(8) 



The purpose of this paper is to give a detailed account of derivation of the revised value 
of the eighth-order coefficient of A\ reported recently [l3] 



-1.914 4 (35) . 



(9) 



Making use o 
divergences Q, 



our automating algorithms in handling ultraviolet (UV) and infrared (IR) 
15] ]. we are now able to generate the eighth-order FORTRAN codes very 
easily and swiftly. However, numerical evaluation of these codes is still nontrivial and 
requires a huge computational resource. Thus far the "new" calculation has achieved a 
relative uncertainty of about 3 % . Although this is still more than an order of magnitude 
less accurate than that of Ref. [l6], it is good enough for the purpose of checking the old 
calculation. 

Comparison of the "new" numerical result with the old one has revealed an inconsistency 
in the treatment of the IR divergence in the latter. With this error of the old calculation 

(8) 

corrected, we now have two independent determinations of A\ . Of course, precise evaluation 
of all terms of "new" A^ by the integration routine VEGAS [ITJ requires an enormous 
amount of computation. Fortunately, as is described in Sec. HVDt the correction term itself 
can be evaluated easily and very precisely. This is why we are able to give the uncertainty 
in Eq. (Q which is essentially identical with that of the previous calculation 16| . 



Besides A\ the known terms of Eq. (CO) are as follows 



Q,Q, 19, 20, 21, 3 



23 



24 



.3 



26|,l27|, |2g,l2S|: 



(2) 



(4) 



(6) 



0.5, 

-0.328 478 965 579 • • • , 
1.181 241 456 587 - ■• , 

4 10) = 0.0 (4.6), 

A^\m e /m^) = 5.197 386 70 (27) x 10" 
Af{m e lm T ) = 1.837 63 (60) x 10" 
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Ai'{m e /m,j) 
Af\m e /m T ) = -6.581 9 (19) x 10~ 8 , 
Af\m e /m^m e /m T ) = 0.190 945 (62) x 10~ 12 



-7.373 941 58 (28) x 10 



-6 



(10) 



Here, A^ \ A± , and A ® are known analytically. A^\ A% and Af^ are known analytically 
as functions of mass ratios so that their uncertainties are due to those of measured lepton 
masses only. Note that A^ is actually unknown and the value listed above is an educated 



guess calculated by the recipe proposed in Ref. [29j to indicate a likely range of the value 
taken by A^. This will soon be replaced by a real number, which is being evaluated 
by FORTRAN codes generated with the help of the automatic code generator GENCODEiV 
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, Jjij]. Until then A^ in Eq. (flOj) is the largest source of theoretical uncertainty. 
In order to obtain the numerical value of the theoretical g — 2, an explicit value of the 
fine structure constant a, which is determined by the physical phenomena other than g—2, 
is required. At present the best values of a available in the literature are from the Cesium 



atom experiments 



30, 



3l| and the Rubidium atom experiment 32 ] 



(11) 
(12) 



a _1 (Cs06) = 137.036 000 00 (110) [8.0ppb] , 
a _1 (Rb06) = 137.035 998 84 (91) [6.7ppb] . 

They lead to the theoretical predictions of a e : 

a e (Cs) = 1 159 652 172.99 (0.10)(0.31)(9.32) x 10~ 12 , 

a e (Rb) = 1 159 652 182.79 (0.10)(0.31)(7.71) x 10" 12 , (13) 

respectively, where the uncertainty 0.10 comes from the eighth-order result ([9]), 0.31 is an 
estimated uncertainty of the tenth-order term, and 9.32 and 7.71 come from the uncertain- 
ties of the input values of the fine structure constants given in Eqs. ( TTT1) and f|T2l) . The 



uncertainty due to the hadronic and weak contributions ([6]) is 0.02 x 10~ 12 . The revised 
theoretical anomaly a e is in closer agreement with the experimental values ([2]) and ([3]) than 
the old value [16 ]. 



Unfortunately, the precision of a e given in Eq. (1131) is not high enough for direct confronta- 
tion between the experimental and theoretical a e 's. This is because the uncertainties in a e 
due to these a's amount to 9.3 x 10~ 12 for a(Cs06) and 7.7 x 10~ 12 for a(Rb06), respectively, 
which are an order of magnitude larger than the experimental uncertainty 0.76 x 10~ 12 and 
the theoretical uncertainty 0.28 x 10~ 12 of a e . 

This implies that, assuming the validity of QED, the electron g—2 is in fact the best source 
of the fine structure constant a, an order of magnitude better than any alternative. Because 
of high precision of the experiments ([2]) and ([3D the fine structure constant a determined 
from a e is rather sensitive to the revision of the theoreticalprediction. Equating the Harvard 



measurements (J^J) or (J3J) , and the theory (jlj), we obtain [71. l33l. \34\ 



a _1 (a e (HV06 = Th07)) = 137.035 999 070 (12) (37) (90) [0.71ppb] , (14) 
a _1 (a e (HV08 = Th07)) = 137.035 999 084 (12) (37) (33) [0.37ppb] , (15) 

(8) 

where the first and second uncertainties come from the numerical uncertainties of A\ and 
A^°\ respectively, and the third in Eq. (JT41) or Eq. (115p comes from the experiment (J2J) or 
(|3J), respectively. 

These values of a -1 are smaller than the old a _1 (a e (HV06 = Th06)) by —6.411 80(73) x 
10~ 7 which is about 4.7 ppb (or about 7 s. d.), but are still in good agreement with a _1 (Rb06) 
of Eq. ( TTTT) and a _1 (Cs06) of Eq. ( fl~2l) . whose uncertainties are about 7 ppb. 

The organization of the paper is as follows. In Sec. II, we briefly overview the "old" 
and "new" approaches to the numerical calculation of the electron g — 2 in QED. In Sec. 
Ill, the diagrams of Group V of the eighth-order term are discussed. We compared the 
results of the "old" and "new" calculations and found an unaccountable difference in the 
results of the diagram M\%. In Sec. IV, the diagram is closely examined instead of 
Mi§. This is because Miq has a similar structure to Mis, but somewhat simpler. We found 
a source of the discrepancy between the "old" and "new" results and the errors in the "old" 
calculation of Mi 6 and Mi 8 are corrected. Sec. V gives the summary of the updated value 
of the eighth-order contribution to the electron g—2. 

Appendix A presents the tests of the automation system GENCODEvVfor the fourth-order 



and sixth-order g — 2's. Appendix B gives our renormalization scheme of the magnetic 
moment amplitude in the "new" approach. Similarly, Appendix C gives the renormalization 
scheme of the renormalization constants in the "new" approach. 



II. OLD VS NEW APPROACH 

(8) 

The purpose of this paper is the presentation of the results of evaluation of A\ by two 
independent methods. Although these methods started from the same Feynman-parametric 

(8) 

representation of A\ , they took different approaches, in particular, in the handling of the 
self-energy subdiagrams and associated infrared (IR) divergences. Furthermore, the "new" 
approach was instrumental in discovering an error in the handling of infrared divergence in 
the old method [16]. After correcting this error, we now have two independent evaluations 
of A\ , enhancing substantially the credibility of the calculation. 



A. Common starting point 

The anomalous magnetic moment a e is given by the static limit of the magnetic form 
factor that is related to the proper vertex part IT. Throughout this paper our attention 
is focused on the q-type diagrams, namely, proper vertex diagrams that have no closed 



lepton loops. In both o 
Takahashi identity 



35 



d and "new" formulations, we use a relation derived from the Ward- 



36] 



A u (p,q) ~ -q 



dA^(p,q) 
dq u 



aE < p) (16) 



between the self-energy part £(p) and the sum of vertex parts A u (p, q) obtained by inserting 
an external vertex in the lepton lines of £ in all possible ways. Here, the momentum of 
the incoming lepton is p — \q and that of the outgoing lepton is p + \q. By means of 
Eq. (fl6l) a set of vertex diagrams are amalgamated into a single self-energy-like diagram, 
which reduces the number of independent integrals substantially. For the eighth-order q- 
type diagrams, the number of Feynman diagrams is reduced from 518 to 74. Taking into 
account the time-reversal symmetry, the number is further reduced from 74 to 47. 

The amplitude of the magnetic moment contribution of a diagram is obtained by apply- 
ing Feynman-Dyson rules of QED in the momentum space. Carrying out the momentum 



6 



integration analytically, we can express the amplitude of 2nth-order diagram Q as an integral 
over the Feynman parameters zf 



M, 



(2n) 



-- (n-1)! (dz) 



1 ( E + C Ei + Ci 
_n - 1 Vt/ 2 ^™" 1 U 3 V n ~ 2 

,'Nq + Zq N l + Z 1 



(17) 



Ipyn JJ3yn-l 

where (dz)g = f^dz^l — ^i 2 *)- The factor (a/7r) ra is omitted for simplicity. 

The quantities E^, Ck, Nk, and Zf, are polynomials of symbols called building blocks 



Bij, Aj, and Cjj 



351 ] . The symbols B^ and C/ are homogeneous polynomials of Feynman 



parameters, related to the flow of loop momenta in the diagram. The symbol Ai is called 
scalar current that is associated with the flow of external momenta. They are functions of 
Bij, U, and Z{. The symbol CV, is given by Zi,Bij and U. The denominator function V is 
defined by 



y = Y j z i -G, G = J2^ 



(18) 



where the summation is over the electron lines only, and the electron mass is chosen to unity 
for simplicity. 



B. Different structure of integrand 

Although the "old" and "new" methods have the common starting point, they have an im- 
portant difference in practice. In the "old" version of the programs, the size of the integrand 
was reduced by taking symmetries of a diagram into account. One type of modifications was 
applied to the integrand by using 8 "junction laws" and 4 "loop laws" satisfied by the scalar 



currents Ai (where Feynman parameters z^ play the role of resistance) [35|, |36(. Another 
type of modification was to reduce the number of integration variables by exploiting the 
fact that in some diagrams the integrand depends only on a particular combination of Feyn- 
man parameters. These resulted in substantial reduction in the size of integrands and the 
amount of computing time required to achieve desired precision. In the "new" version, those 
modifications were not employed at all because they are diagram-specific and not suitable 
for automation. As a result, the size of FORTRAN source code for M i (see Fig. [[]), which 
requires only vertex renormalization, is about 515KB in the "new" version in contrast to 
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316KB of the old version. A more notable difference is seen for the diagram M\i which 
requires IR subtraction. The "new" M12 occupies 630KB while the old Mi 2 occupies only 
21KB. As a consequence, the old and "new" integrals have much different forms so that they 
can be regarded to be independent of each other as far as numerical integration is concerned. 



C. Ultraviolet (UV) divergence 



The amplitude Mg has (logarithmic) UV-divergences in general. Suppose we want to find 
out whether Mg diverges when all loop momenta of a subdiagram S consisting of N$ lines 
and n$ closed loops go to infinity. In the parametric formulation this limit corresponds to 
the vanishing of U when all Zi for % e S vanish simultaneously. To find the criterion for the 



UV divergence from S, consider the part of t 
X^ies z i — e - m the limit e — > one finds 36 



re integration domain where Zi for % e S satisfy 



V = 0{1), U = 0{e ns ) 

0{e ns ' 1 ) for i,jeS, 
0(e ns ) otherwise. 

The UV-divergent part can be identified by the following procedure called .fC-operation: 



Bij 



(19) 



(a) In the limit (1191) keep only terms with the lowest power of e in U, By, and Aj. In this 
limit U factorizes as UsUg/s where Q jS is obtained from Q by shrinking S to a point in 
Q. Bij factorizes similarly. V is reduced to Vg/s, where Vg/s is the V function defined 
on Q/S. 

(b) Replace Vg/ S by Vg/ S + V s . 

(c) Rewrite the integrand of Mg in terms of parametric functions redefined in (pjj and (jbj), 
and drop all terms except those with the largest number of contractions |35[] within S. 
The result is denoted by KgMg, in which Kg stands for an operator acting on Mg. 

By construction, K$Mg has the same UV divergence as Mg in the same integration domain. 
Therefore it can be used as a pointwise subtraction term in the subtractive renormalization. 

An important feature of ii'-operation is that the resulting integral can be factorized exactly 
into a product or a sum of products of lower-order quantities that consists of a leading 



8 



UV-divergent part of the renormalization constant and the magnetic moment part. The 
-ff-operation associated with a UV- divergent subdiagram S produces, when S is of vertex 



type, the subtraction term of the form [36j, |37] 



K s Mg = LfMg/s , (20) 

where L^ v is the leading UV-divergent part of the vertex renormalization constant Ls and 
Mg/ S is the magnetic moment part of the reduced diagram Q/S. When S is a self-energy 
subdiagram, the i^-operation yields 



36 



K s M g = Sm^Mg/s^) + BfM g/m , (21) 

where 5m^ v is the leading UV-divergent part of the mass-renormalization constant 5ms, 
is the leading UV-divergent part of the wave-function renormalization constant Bs, and 
the reduced diagram Q/[S, i] is obtained from Q by removing S and a lepton line i adjacent 
to S. 

The whole UV-divergent structure of the amplitude Mg can be recognized by Zimmer- 

L 

mann's forest formula 14]. A forest is a set of UV-divergent subdiagrams in which any 
pair of subdiagrams is either disjoint (they do not share lines or vertices) or inclusive (one 
subdiagram is a subgraph of the other subdiagram). Each subtraction term corresponds 
to a forest. In our formulation, the subtraction term is obtained by successive application 
of i^-operations for every element of the forest. The UV- finite amplitude M_g created by 
./^-operation is thus expressed in the form 

M g = Mg + J2U ("Ks) Mg , (22) 

/ 56/ 

where the summation is taken over the normal forests of the diagram Q that do not include 
Q itself as an element. 

N. B. In both old and "new" approaches UV divergence is treated by the same K-operation. 



D. Infrared (IR) divergence 

A diagram may have an IR singularity when some of the internal photon momenta vanish. 
In order that this singularity becomes actually divergent, however, it must be enhanced by 
vanishing of denominators of two or more electron propagators (called enhancers) due to 

9 



kinematical constraints. Such a situation occurs in diagrams that have self-energy-like sub- 

re vanishing of l^-function of the denominators 

3,0 



diagrams. In Eq. ((TTJ) this corresponds to t 
in the integration domain characterized by 



0(5) if i is an electron line in 7Z, 

(9(1) if i is a photon line in 7Z, (23) 

0(e), e~5 2 , itieS, 
where 7Z = Q/S. 

This enables us to obtain a simple IR power-counting rule for identifying IR divergent 
terms. When there are two enhancers, the amplitude shows a logarithmic IR divergence. We 
can identify and construct the corresponding subtraction term by the following procedure 



called 3-operation [36l. l37j|: 



(a) In the limit (1231 keep only terms with lowest power of e and 5 in 11,73^, Ai. The 
numerator then factorizes to the product 

F -> F n F s , (24) 

where F-ji is a numerical factor obtained by replacing all scalar currents Aj in the diagram 
7Z by one. 

(b) Make the following replacements: 

U^U s U n) V^Vs + Vk, F - F [L n ]F s , (25) 

where Fo[Ltj,] is the no-contraction part of the vertex renormalization constant defined 
in 71. The difference between Fo[L^] and F-ji causes a finite difference of the integration. 

(c) Rewrite the integrand of Mg in terms of redefined parametric functions, keeping only 
the IR-divergent terms. 

In the "old" method all logarithmic IR divergences have been subtracted by means of the 
3-operation. However, the case involving linear IR divergence, which has three enhancers, 
was handled by an ad hoc manner instead of a systematic approach. Actually, the cause 
of linear IR-divergence is easy to identify. It is caused by our treatment of self-energy 
subdiagram by means of i^-operation which subtracts only the UV-divergent part of the 

10 



self-mass. The unsubtracted part of self-mass keeps the number of enhancers unchanged, 
except in second-order case where the /T-operation subtracts the self-mass term completely. 

In the "new" approach, a systematic method is developed to handle the linear IR di- 
vergence. To remove the finite remnant of self-mass term completely, an /2-subtraction 
operation 



15| is newly introduced. After the i?-subtraction operation is carried out, which 
decreases the number of enhancers to two, only logarithmic IR divergences remain, which 
can be handled by the /-subtraction operation similar to, but different in detail from, the 
3-operation of the "old" method. 

For a formal treatment, we introduce two operators for these subtractions. The R- 
subtraction operator M.$ acts as 

M s M g = 5m s i M g/s{i , ) , (26) 

where 5m^ is the residual part of the mass renormalization constant defined by 

5m^ = 5m s - 5m^ Y + (~ K s') (27) 

/ s'ef 

in which the leading UV-divergent part 5m^ v and the subdivergent parts associated with 
the forests Ils'e/ ( — ^5') 5ms are subtracted away, where 5m = 5m — <5m uv . 
The /-subtraction operator Is acts on the UV-renormalized amplitude M_g as 

hK g = L* /s Ms , (28) 
where Lg, s is the residual part of the vertex renormalization constant defined by 

L^ /s = L G/S - L^ g J s + (" K 5') l g/s (29) 

/ S'ef 

in which the leading UV-divergent part LgJ s and the subdivergent parts associated with the 
forests ris'e/ ( — ^s') Lg/s are subtracted away, where L = L — L uv . 

N. B. The IR power counting rule identifies only IR-divergent terms. It does not specify 
how to handle IR-finite term. The "new" I-subtraction operation handles the IR-finite terms 
differently from the "old" ^-operation. The I-subtraction operation needs not deal with the 
IR divergence associated with a vertex subdiagram of the self- energy- like diagram, while the 
3 -operation directly acts on the vertex subdiagram. 

The whole set of IR subtraction terms can be obtained by the combinations of these two 
operations, both of which belong to annotated forests [3] . An annotated forest is a set of 
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self-energy-like subdiagrams, to each element of which the distinct operation of /-subtraction 
or .^-subtraction is assigned. The IR-subtraction term associated with an annotated forest 
is constructed by successively applying operators I or K, and takes the form 

(-!&)...(-%.)... Me (30) 

where the annotated forest / consists of the subdiagrams Si, ... and Sj, .... 

N. B. The IR divergence is treated differently in the old and "new" approaches. This 
difference plays an important part in ensuring the independence of two calculations. 



E. Residual renormalization 



Because of difference in the handling of IR divergences in the "old" and "new" methods, 
we obtain different forms of residual renormalization. Since the "old" residual renormaliza- 
tion is described in Refs. 



36 



37l |. let us consider here only the "new" residual renormaliza- 



tion. 



In the "new" approach the UV- and IR-finite amplitude has the form 



AMg = Mg + J2 Y[(-K S )Mg + 

/ 56/ J 



Q ■ 



(31) 



where M g is the UV-fmite quantity defined by E q. f[2"2"j) . AMg can be readily turned into 



14 



151 ] and is to be evaluated by numerical 



a numerical integration code by gencodetV 
means. 

This procedure is different from the standard on-shell renormalization which is defined 
by the on-shell quantities. The difference between the on-shell quantities and the quantity 
evaluated by Eq. (l3Tj) must be compensated by products of known lower-order quantities. 
We call this step the residual renormalization. See Appendix B for details. 



III. EIGHTH-ORDER TERMS 

The eighth-order term A\ receives contributions from 891 Feynman diagrams. The 
373 of them have closed lepton loops and had been evaluated by two or more independent 
methods [16J. The remaining 518 diagrams of g-type form one gauge-invariant set (Group V). 
In our approach they are represented by 47 independent diagrams shown in FigJT] by using 

12 



/OO, /toQ /^sTX //T^N /t<7^\\ /yT7^\ 

MOl M02 M03 M04 M05 M06 M07 



M15 



M22 



M29 



M36 



_ 



M43 



M09 



MIO 



M16 




M30 



M37 



M17 



M18 



M19 



^ 

M23 M24 M25 M26 



M31 



M32 



M33 



M38 



M39 



M40 



/f^\ 

M44 M45 M46 




M47 



M14 



M20 



M21 



M27 



M28 



M34 



M35 




M41 



M42 



FIG. 1: Eighth-order Group V diagrams represented by 47 self-energy-like diagrams M01-M47. 

the relation derived from Ward-Takahashi identity and the time-reversal symmetry. Thus 
far, there is only one complete evaluation of the eighth-order term, which was performed 
by numerical means 16]. Some of these diagrams have linear IR divergence, which was 
treated by an ad hoc subtraction method. In contrast GENCODEiV is capable of dealing with 



such hard IR divergence in a systematic fashion [15[ . The application of GENCODEiV to the 
calculation of the eighth-order (/-type diagrams provides us the opportunity not only to test 
if it works properly, but also to check the previous result. 

Even in the eighth-order case GENCODEiV creates FORTRAN programs very rapidly. The 
entire 47 program sets are generated in less than ten minutes on hp's Alpha. The numerical 
evaluation is, however, quite non-trivial and requires a huge computational resource. For 
the preliminary evaluation we have used 64 to 256 Xeon CPU's per diagram and run the 
programs over a few months. To our surprise it uncovered an inconsistency in the treatment 
of IR subtraction terms in the old calculation. In Sees. II I II and IIVI we describe how this 
inconsistency was uncovered by a detailed comparison of the old code and the code generated 
by GENCODEiV. 



13 




FIG. 2: Self-energy- like diagrams of M%q and M\s- Feynman parameters assigned to the electron 
lines are z\ ~ zj and those to the photon lines are z a ~ z&. 

A. IR treatments of the eighth-order diagrams 

The treatments of IR subtraction terms are different in gencodeA^ and the "old" ap- 
proach. The difference of IR subtraction terms leads to the difference of the finite part 
of the amplitude AM; (i = 01, ■ • • , 47). The difference AM° ld - AMf w between the old 
amplitude AM° ld and the "new" one AM™ is finite and can be expressed analytically in 
terms of finite lower-order quantities. We will see if this difference is numerically reproduced 
by substituting the numerical values calculated separately for these lower-order quantities. 
If the numerical discrepancy is found, there must be something wrong in either the "old" or 
"new" calculation. This is what we tried to find out. 

We noted in Sec. Ill CI that the subtraction of UV divergences is achieved by the same 
A"-operation in both GENCODEiV and "old" calculation. Therefore, the difference between 
AM° ld and AMf ew , if it exists, comes exclusively from the difference of IR subtraction 
procedures between the "new" and "old" calculations. To examine this difference more 
closely let us begin by considering relatively simple diagrams which contain only one second- 
order self-energy subdiagram. Diagrams belonging to this class are M02, M03, M09, M13, 
M14, Mi 5 , M23, M 24 , M27, M43, M44. As an example let us consider the diagram M i4 . 

The IR divergence occurs in M14 from the second-order self-energy subdiagram which 
consists of an electron line "2" and a photon line "b" in Fig. [3]^a). In the W-T summed 
diagram this subdiagram plays dual roles. One part of this subdiagram behaves as a genuine 
self-mass term and the associated UV singularity is removed completely by the K 2 -operation. 
Another part works as the second-order magnetic moment M 2 , and the residual diagram 
surrounding M 2 behaves like a sixth-order vertex diagram of Fig. E](b), which is IR- 
divergent. 
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FIG. 3: (a) Self-energy-like diagram of M\ 4 . (b) Vertex diagrams 6g5 and 41 from which the vertex 
renormalization constants Lq 9 ^ and L 4 i, respectively, are derived. 

In the "old" approach, the finite contribution AM U was defined by 

AM°4 d = M 14 + l[(-^s)M 14 - (J 6ff5 + I 2 AL 4l )M 2 . (32) 
/ sef 

Here the second term on the right-hand side is the sum of UV subtraction terms given by 
the if-operation. The last two terms beginning with the letter "J" are the IR subtraction 
terms generated by the 3-operations J134567 and ^13(1 — J134567) in the "old" approach. They 
arise from the "magnetic-moment part" of the self-energy-like subdiagram mentioned above. 
Note that they are exactly identical with the IR-divergent parts of L 6g5 : 

+ I 2 AL 4l + AL 6g5 - Y[(-^s)L 6g5 + L^J . (33) 

/ sef 

Here the sum appearing on the right-hand side denotes all the UV sub divergences contained 
in L eg5 (whose explicit form is — L\ v L 4c — L^ V L 2 + (L% V ) 2 L 2 ), and the last term is the 
overall UV divergence of L eg5 . 

In the "new" (or GENCODEiV) approach, we introduce a term Lf g5 defined by 

L* 5 = L 695 - LgJ + J2 II(-^) Z 6 9 5 (34) 

/ sef 

and AMJ 1 ™ by 

AM° 4 CW = M 14 - I 2 M 14 , (35) 

where I2 is an /-subtraction operation associated with the self-energy-like subdiagram S = 
{2, b}, and yields 

h M14 = L* g5 M 2 . (36) 
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From Eqs. (l33|) and (|34|) we obtain 

-^695 = -^6g5 + hAL^ + AL Gg5 . (37) 

This is UV-finite and consists of not only IR divergent terms but also a completely finite 
term AL &gb . It follows that AMf 4 d and AM™ differ by 

AM* - AM™ = AL 6g5 M 2 . (38) 

Since AL 6g5 is UV- and IR-finite, it can be computed without encountering with UV or IR 
divergence. This is true for every AM° ld — AM™ as it originates from the choice of finite 
pieces that accompany the singular terms. The choice adopted in AM"| W turned out to be 
preferred since it leads to a simpler formula and can be readily extended to other cases. 

All eleven diagrams listed above can be analyzed in the same manner. The diagrams M04, 
Mn, M12, M i7 , M 2 g, M 30 , which contain two or three second-order self-energy- like subdia- 
grams, are slightly more complicated, but can be treated in a similar manner. Evaluation 
of the diagrams with one self-energy subdiagram of fourth- or sixth-order such as Mos, M10, 
M26, M38, M40, M41 is more complicated and needs the residual self- mass renormalization, 
the /^-subtraction, as well as the /-subtraction. But they do not present particular diffi- 
culty as far as IR subtraction is concerned. (See Appendix B for more information on these 
diagrams.) 

The diagrams M 2 s, M 42 , M 45 , M 46 , M 47 are even more complicated due to nested struc- 
ture, but they can also be handled by slight extensions. (See Appendix B.) 

The most difficult of the eighth-order q-type diagrams are those containing one second- 
order self-energy- like subdiagram and one fourth-order self-energy-like subdiagram, namely, 
Mi6 and Mig of Fig. [2j The difficulty originates from the fact that these diagrams have 
linear IR divergence. Detailed analysis of these diagrams is deferred to Sec. HVl 

B. Numerical result of eighth-order calculation 

We present the results of our numerical study for AM° ld — AM™ in Tables H] HT1 and IITT1 
In these tables, the second columns list the analytic expression of AM° ld — AM™ in terms 
of finite pieces of lower-order renormalization constants and magnetic moment amplitudes 
multiplied by the multiplicity, which is 1 for the time-reversal-symmetric diagram and 2 
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TABLE I: Comparison of the numerical calculation of M01-M15 of the eighth-order Group V 
diagrams. The second column shows the analytic expression for AM° ld - AM° ew for each dia gram 
Mi in terms of lower-order finite quantities multiplied by the multiplicity. The value A in the third 
column is obtained by substituting the values of lower-order renormalization constants, such as 
AM 4a , AL 4s and AM^jp), for the corresponding expression in the second column. See Appendices 
B] and [C] for the detail. In contrast, the value for AM° ld — AM" ew in the fourth column, denoted 
by value B, is obtained by taking the direct difference between the value of AM° ld quoted from 
:lef. Ilal . and the one of AM 4 new calculated via GENCODEiV in the "new" IR subtraction procedure 
LJ, [l5|. The fifth column lists up the differences A — B. If the whole calculation is done correctly, 
A — B must vanish within the numerical uncertainty. In evaluating AM ncw the double precision 
is used for the diagrams without a self-energy subdiagram, while the quadruple precision is used 
for the remainder. 



Diagram 


difference 


value A 


value B 


A - B 


M i 








-0.0129(47) 


0.0129(47) 


M 02 


2AL 6/ iM 2 


-0.0063(2) 


0.0060(110) 


-0.0124(110) 


M 03 


AL 6/3 M 2 


-0.1133(1) 


-0.1055(100) 


-0.0078(100) 


M 04 


2{AL 6dl + AL 6d3 )M 2 


0.3350(2) 


0.3408(175) 


-0.0058(175) 


M 05 








0.0020(28) 


-0.0020(28) 


M 06 








-0.0223(61) 


0.0223(61) 


M 07 








-0.0102(40) 


0.0102(40) 


M 08 


2(A<5m 4 aAM 4a(1 * ) + AL 4c AM 4a ) 


-2.1809(7) 


-2.1790(121) 


-0.0019(121) 


M 09 


2AL 6/2 M 2 


0.0806(1) 


0.0894(109) 


-0.0088(109) 


M w 


2(A<5m 46 AM 4a(1 *) + AL 6d2 M 2 + AL 4c AM 4fe ) 


15.8898(49) 


15.8795(147) 


0.0103(155) 


Mn 


2AL M5 M 2 


0.6949(2) 


0.6827(112) 


0.0122(112) 


Mi 2 


(2AL 6ai + AL 6a3 )M 2 


1.2842(0) 


1.2875(74) 


-0.0034(74) 


M13 


2AL 6hl M 2 


-0.4211(2) 


-0.4238(48) 


0.0027(48) 


M M 


2AL 6g5 M 2 


0.0892(2) 


0.0960(95) 


-0.0068(95) 


Mis 


2AL 6gl M 2 


0.0883(2) 


0.0893(71) 


-0.0009(71) 
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TABLE II: Comparison of the numerical calculations of Miq-M^q of the eighth-order Group V 
diagrams. 



Diagram difference value A value B A — B 



71 /T 


2(Adm4 AM4 6 ( 1 *) + AL 6d M 2 + AL 4s AM 4a ) 


— 2.b042(b) 


— 2.631b(235) 


0.0274(235) 


71 /T 


2(AL 6e i + AL 6(j4 )iW 2 


— z.lzUl(z) 


imn/ion\ 

—2.1010(189) 


—0.0173(189) 


7\ /T 

Mis 


Of A A7V/f iAT A7V/f 


16.9686(39) 


17.1897(20b) 


—0.2207(210) 




i/AT 1 AT s \7l/fl 

+ (AL 6b i + AL 6a2 )M 2 | 








Mig 


u 





U.UUU2(o) 


— (J.U(J(J2(3j 




U 


u 


n nm nVi t\ 
U.UUlU^l ( J 


n nm nV 1 t\ 

— u.uuiu^i r ) 


71 /T 

Mix 


u 


U 


0.0003(3) 


—0.0003(3) 


71 /T 


u 


U 


—0.0090(25) 


0.0090(25) 


71 /T 


A r 71 /f 


0.0501(2) 


0.0438(59) 


0.00b4(59) 




2AL fi „9M 9 


0.0789(2) 


0.0945(61) 


-0.0155(61) 


M 25 








-0.0031(20) 


0.0031(20) 


M 26 


A(5m 6/ (M 2 * -M 2 *[I]) 


2.5119(3) 


2.5369(95) 


-0.0250(95) 


M 27 


2AL 6s4 M 2 


-0.0630(1) 


-0.0459(90) 


-0.0171(90) 


M 28 


2{A8m ed (M 2 * - M 2 *[I}) + AL 6c2 M 2 } 


-7.5332(5) 


-7.5307(153) 


-0.0025(153) 


M 29 


2AL 6e2 M 2 


-0.2857(2) 


-0.2809(109) 


-0.0048(109) 


M 30 


A5m 6a (M2* - M 2 * [I]) + 2AL 6b2 M 2 


0.2763(6) 


0.2675(153) 


0.0088(153) 



otherwise. Each value in the third columns, called "value A", is obtained by substituting 
the values of these renormalization constants, etc., listed in Table IV, for the corresponding 
expression in the second columns. 

In contrast to value A, each value in the fourth columns, called "value B", is obtained 
by taking the difference between the numerical value AM° ld quoted from the literature [3] 
and the one AMf ew newly calculated via gencode A according to the "new" IR subtraction 

n 

procedure [15(. The fifth columns list up the difference of value A and value B for each i, 
denoted by A — B. It must be zero within numerical precision if the whole calculation has 
been done correctly. If value A and value B are different, there are two possible sources. 
One possibility is that the program used for a numerical calculation has a bug. It means 
that either AM* or AMf ew is wrong, or both are wrong. The other possibility is that we 
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TABLE III: Comparison of the numerical calculations of M31-M47 of the eighth-order Group V 
diagrams. 



Diagram 


difference 


value j4 


value B 


A-B 


M31 








0.0007(5) 


-0.0007(5) 


M 32 








-0.0024(10) 


0.0024(10) 


M 33 








0.0001(3) 


-0.0001(3) 


M 34 








-0.0010(13) 


0.0010(13) 


M35 








0.0001(13) 


-0.0001(13) 


M 36 








-0.0027(22) 


0.0027(22) 


M 37 








0.0004(5) 


-0.0004(5) 


M 38 


A6m 6h (M 2 * -M 2 *[I] ) 


-0.9088(3) 


-0.9112(40) 


0.0024(40) 


M 39 








-0.0031(18) 


0.0031(18) 


M 40 


2A<5m 69 (M 2 * -M 2 *[I] ) 


3.8281(3) 


3.8326(71) 


-0.0045(71) 


M41 


A5m 4a (AM 4a(2 *)) + AL 4x AM 4a 


0.9809(3) 


0.9713(83) 


0.0096(83) 


M 42 


A5m 6c {M 2 * - M 2 *[I]) + AL 4i AM 4a 


-7.0218(4) 


-7.0202(114) 


-0.0016(114) 




+A<5m 4 a{AM 46(2 * ) - A5m 2 *(M 2 * - M 2 *[I])} 








M 43 


AL 6h3 M 2 


0.4724(1) 


0.4703(42) 


0.0022(42) 


M 44 


2AL 693 M 2 


-0.0748(1) 


-0.0499(69) 


-0.0250(69) 


M 45 


A5m 6e (M 2 * - M 2 * [I]) + AL 6c3 M 2 


-0.0523(3) 


-0.0498(90) 


-0.0025(90) 


M 46 


A5m 46 AM 4a(2 *) + AL 6e3 M 2 + AL 4x AM 4fe 


-7.9339(22) 


-7.9232(86) 


-0.0107(89) 


M 47 


A5m 66 (M 2 * - M 2 * [/]) + AL 663 M 2 + AL 4 jAM 46 


10.5872(15) 


10.5864(102) 


0.0008(103) 




+A<5m 46 {AM 4b(2 * ) - A5m 2 *(M 2 * - M 2 *[I])} 









incorrectly identified the analytic difference between the "old" and "new" methods. 

For a diagram Mj without any self-energy-like subdiagrams, the analytic expression of 
AM° ld — AM" ew is trivially zero, as it does not have IR divergence. We can see from the 
corresponding values B in Tables [H [III and II III that this is confirmed within the numerical 
precision employed. 

The diagrams containing self-energy-like subdiagrams suffer from IR divergence. Tables[H 
HT1 and II I II show that "old" and "new" calculations are in good agreement for most of these 
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TABLE IV: Finite renormalization constants used in Table HI [UJ an d IIII1 Sixth-order vertex 
renormalization constants are shown in this table. Their validity is checked by comparing the sum 



Xlbd_ = Ei=i AL 6« + iABfc + 2A5m 6x , 



x 



h to the previous Xlbd values listed in 



Ref. [jj. 


AL« i 


539604( 45) 


ALr n 


—0 16721K 81) 


ALr q 


1 489159( 98) 




— 1 479745(109') 


A-Zyg{,2 


582944(106) 


A Lrkq 


—0 016344( 73) 


ALr i 


— D 219365 ( 98") 


ALr i 


071504( 87) 


ALr q 


—0 552261(107) 


A Lrji 


n 834949 ( 96") 




—0 090796( 92) 


A-t/6d3 


—0 499995( 97) 


A TjRJA 


— 1 378190(109"! 




694916(101) 






A J 

^■^bel 


—0.741904(144) 


A J 


—0.285670(108) 


A J 

^-^6e3 


r\ -i /I 1 '70'7/l on) 

—0.141787(122) 




-0.006322(114) 




0.080648( 97) 


AL6/3 


-0.226693(106) 


AL 6gl 


0.088204( 70) 


AL 6s2 


0.078922(103) 


A-^693 


-0.074834( 92) 


AL 69 4 


-0.062995( 85) 


AL 6g5 


0.089213( 69) 








-0.421132(108) 


AL6/i2 


0.050140(108) 




0.944887(116) 


A5m 6a 


-0.15331(26) 


A£m 6b 


1.83795(19) 


AJm 6c 


-3.05047(17) 


A(5m 6(i 


-1.90117(11) 


A5m6e 


0.11193(13) 


A^mgj 


1.25594(10) 


A<5m 69 


0.95702(6) 




-0.45441(5) 







diagrams. However, a large discrepancy —0.221 (21) is found for the diagram M\%. Though 
no detectable discrepancy is found for Miq, it has a structure similar to Mi§ and is somewhat 
simpler to analyze. In Section [IV] we thus look for the origin of such a discrepancy through 
a detailed investigation of Mi 6 . 

IV. DETAILED EXAMINATION OF M 16 



In the "old" approach the finite contribution AM 16 was given by 36, 



f sef 



- IbdM 2 - -Jb C M 2 - h s AM Aa - A5m ia hb{i*) + h*^m Aa M 2 



(39) 
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TABLE V: Finite renormalization constants used in Table |IJ |TT] and IIII1 Fourth-order and second- 
order quantities are given here. 



ALa 


n 003387('16 s i 


ALa 


—0 481834(54'! 


ALa 


407633(20) 


A Lai 


124796(67) 


A R. 


— u.uoyoi i^io j 


AR„ 


^Q79R^/' 1 t\\ 
— U.oy i ZooylQ j 


A<5m 4a 




/\ A T K 1 . i 

L\(J 1 1 LAJ) 


Z.ZU i / I I J 


A M, 


u.zioooy ^oy j 


A M,, 


— u.ioi ozo^oy j 


AM 4fl(11) 


3.6192(31) 


AM 4a(2 *) 


-3.6003(19) 


AM 4b(1 *) 


4.2486(15) 


AM 46(2 *) 


1.6432(15) 


AM 2 


0.5 


AM 2 * 


1 


AM 2 *[I] 


-1 


A5m 2 * 


-0.75 


AB 2 


0.75 


AB 2 * [I] 


-0.5 


AL 2 * 


-0.75 


AB 2 * 


1.5 



while the "new" version is given by 

AM™ = Mia + X] n(- K 5)M 16 

/ 56/ 

- L* cl M 2 - LjAM 4o - A5m 4a M 46(1 , ) + L\\A5m 4a M 2 , (40) 

where Lf* = I 2 *. Note that "2*" denotes the second-order diagram with a two-point vertex 
inserted into the internal lepton line. "46(1*)" denotes the diagram obtained from the 
fourth-order diagram 46 by inserting a two-point vertex into the lepton line 1 . 

A. Unrenormalized amplitude and UV subtraction terms of Miq 

We began our examination by comparing the unrenormalized amplitude M\§ and its UV 
subtraction terms in the "old" and "new" programs. For this purpose we used the "spot- 
check" method, by which the values of "old" and "new" integrands are compared at the 
same set of numerical values of integration variables. The integrand of Miq is defined in the 
Feynman parameter space that spans a hyperplane in 11-dimensional space satisfying 

z 1 + z 2 H h z 7 + z a + z b + z c + z d = 1 . (41) 
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In the "new" version, this hyperplane is mapped onto a unit 10-dimensional hypercube. On 
the other hand, in the "old" version the integration space is mapped onto an 8-dimensional 
hypercube, since the integrand depends only on the combination of Feynman parameters 
^137 = Z\ + z 3 + z 7 . To carry out the spot check, we must use the same mappings, so we 
changed the mapping of the "new" integrand to the "old" one defining z\ = z 3 = zj = 
(1/3)^137. In practice the set of input parameters is chosen from the neighborhood of the 
singular point of interest where numerical disagreement is likely to be magnified. But points 
too close to the singular point are avoided, where the noise due to round-off error obscures 
the meaningful information. The "old" integrals and "new" integrals of the unrenormalized 
term and UV subtraction terms should be algebraically equivalent but have different forms 
because of extensive simplification of the "old" integrands by means of various relations 
among scalar currents. The "spot-check" comparison of "old" and "new" unrenormalized 
and UV integrands proves unambiguously that they have nevertheless the same values within 
the precision of numerical evaluation. (Typically more than 10 digits in 14 digits precision.) 



B. IR subtraction terms of M\§ 

The "spot-check" method, however, is not directly applicable for comparison of "old" and 
"new" IR subtraction terms, because they are algebraically different by construction. For 
this purpose we need to understand precisely the analytic structure of the IR subtraction 
terms in both "old" and "new" methods. Thus we follow an alternative approach by which 
we can identify how they differ from each other in the analytic form. 

In the "old" method, the IR singularities of Miq are isolated by the D^-operations, where 
1Z = {1, 3, 7, a}, {1, 2, 3, 7, a, b}, or {1, 3, 4, 5, 6, 7, a, c, d}. Thus, the formal expression of the 



IR-free contribution of M 16 is given by [371 ] 



AM 16 = (1 - J 137 )(l - 3l34567)(l - 3i2 37 )M 16 , (42) 

where M 16 is the UV-finite amplitude obtained by the iT-operations. The product 
^134567^1237 gives no contribution, since they overlap each other and cannot take these IR 
limits simultaneously. 



Following the "old" prescription in Ref . 37| , individual IR subtraction terms of Eq. (142]) 
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can be written as follows: 

" 3l34567Ml6 = -hclM 2 , (43) 

- 3i23 7 M 16 = -J 4c AM 4a - M 46(1 *) [/] A5m 4a , (44) 

- 3x37(1 " 3l34567)(l - ^1237)Ml6 = +/ 2 *A5m 4a M 2 . (45) 

/6ci and J 4c are the no-contraction terms of the vertex renormalization constant Lq c \ and 
L 4c , respectively. Mm\*\ is the magnetic moment amplitude, which is obtained from the 
fourth-order diagram with the two-point vertex inserted into the fermion line 1. Its 
argument [I] implies that the numerator of the no-contraction term of M 4 6(i*) is replaced by 
that of the vertex renormalization constant L 4s , while discarding the contraction terms. 



If there were only logarithmic IR divergence, AM 16 defined in ( T4"2"j) would be IR-finite, 
but it is not. The problem here is that I§ c \ and M 4 &(i*)[I] in Eqs. ( 143]) and (14"4"I) have linear 
IR divergence. The 3-operation prescription is constructed so that it only deals with the 
leading IR singularity. Of all eighth-order g-type diagrams, M 16 and M 18 have the linear IR 
divergence. Since these are the only cases, we chose to deal with their next-to-leading-order 
IR divergences by an ad hoc method rather than constructing a general rule. 

In the 3i3 4 56 7 limit of Eq. f j4"3l) , the diagram Mi 6 decouples into the vertex diagram L 6c i 
which consists of lepton lines 1, 3, 4, 5, 6, 7 and photon lines a, c, d and the magnetic moment 
part M 2 which consists of lepton line 2 and photon line b. All IR singularities originate from 
the vertex diagram L 6cl . The no-contraction term L 6cl [F ], namely J 6c i, includes the leading 
linear IR singularity as well as the next-to-leading logarithmic singularity. 

The logarithmic IR singularity also arises from the 3i37-limit of the one-contraction term 
L§ci[Fi\. To deal with this, we constructed the quantity junrenorm. j n ^he numerator 

is the 3i37 limit of L 6cl [F!], but the denominator V and U is the same as L ecl . The UV 
divergences of j™ renorm - are removed by the 1K 456 , K 56 , K 45 -operations: 

\ J 6c = (1 - K 456 )(l - K 45 )(l - K 56 ) f~J {dz) G/S j^j , 

fi = -16[fi 45 (2 - A 6 ) + 25 46 (1 - 2A 5 ) + B Fj6 (2 - A 4 )] , (46) 

where S = {2, b}. 

Next we consider the 3i237-limit of Eq. (14*41 . In this limit, the self-energy- like subdiagram 
consisting of lepton lines 4, 5, 6 and photon lines c, d plays dual roles. When this fourth- 
order self-energy- like subdiagram behaves as a magnetic moment M 4a , the residual diagram 
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resembles a vertex diagram L is . Its singularity is logarithmic, so Ji237-operation properly 
works for this part. 

The problem arises when the self-energy- like subdiagram acts as the self-mass 5m 4a , which 
is the second one of its dual roles. The residual diagram is the magnetic moment amplitude 
with one two-point vertex inserted, namely, Mm\*)- The power counting shows that it has 
a linear IR singularity. Thus, the J-operation is not enough to remove the IR singularity of 
this term. Mmi*\[I] is not sufficient to remove all IR singularities arising in the J1237 limit. 
The IR structure of M^i*) is more closely scrutinized in the next subsection, where 1^(1*) is 
constructed to include both linear and logarithmic IR singularities of the magnetic moment 
amplitude M^n*)- (A similar subtraction method works also for M 18 .) 



Taking these considerations into account, we rep 
the "old" method listed in Eqs. (g3]) and (jHJ) with 



ace the IR subtraction terms of M 16 in 
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- 3'l3456 7 Ml6 = '(he! + \-hc)M 2 (47) 

- 3I237M16 = -/4 C AM 4a - I 4bil ± } A5m 4a (48) 

which are more convenient for comparison with the "new" approach. Note that Eq.( l4"5l) is 
unchanged. 

Now, let us look at the "new" approach. All IR singularities, both linear and logarithmic, 
are subtracted by using the general rule applicable to any order of the perturbation theory. 
The R- and /-subtractions, and their combinations determine the IR subtraction terms of 
Mi6 as follows: 

-I2M16 = M 2 > 

-l4 56 M 16 = -Ll M^, 
-KU 56 M 16 = -6ml M 4b{1 , } , 

+I 2 M 456 M 16 = +L* 5mf a M 2 . (49) 
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By definition given in Eqs. fl27|) and (|29|) . the residual quantities are explicitly given by 

Lf cl = (1 - K 456 )(l - K 45 )(l - K 56 )(L 6cl - L uv ) 

= L 6cl - L™ - (B™Z' 2 + 5m™L 2 *) - 2L™L is + 2L™(5m 2 L 2 * + B™1' 2 ), 
Ll = (1 - K 2 )(L 4s - L U J) = L 4s - L U J - (^ y Z; + 5m 2 L 2 .), 
M 4a = (1 - K 45 )(l - K 56 )M 4a = M 4a - 2L 2 ^ y M 2 , 
Smf a = Sm Aa - 5m™ 
M«(i*) = (! - K 2 )M 46(1 *) = M 46 (!*) - (5 2 ^M 2 * + 5m 2 M 2 „) , (50) 

where L = L — L . (See below Eq. ([2"T]) .) In terms of the "old" expression of the 
unrenormalized amplitude and renormalization constants, the residual quantities are related 
to the IR divergent and finite pieces of the "old" method by the following relations: 

-^6cl = ^6cl + 2^6c + AL 6cl , 
L^s = I 4s + AL^, 

M 4a = AM 4a , 
Smf a = ASm ia , 

M46(l*) = J 46(l*) + AM 4fe(1 *) . (51) 

We are now ready to compare the IR subtraction terms of "old" and "new" method side 
by side: 

old new 

(a) -(/ 6cl + \ J 6c )M 2 -(/ 6cl + \ J 6c + AL 6cl )M 2 

(b) -/ 4s AM 4fl -(/ 4s + AL 4s )AM 4a (52) 

(c) -J 4b(1 *)A5m 4a + AM 4fe(1 *))A(5m 4a 

(d) +/ 2 *A<5m 4a M 2 +/ 2 *A5m 4a M 2 . 

Actually, instead of examining the IR subtraction terms of the "old" method themselves we 
reconstructed them from the "new" programs by dropping finite terms (eg. AL 4s ) from the 
"Residual" term (eg. Lf s ), and compared them with the terms in the "old" programs by the 
spot-check method. To obtain Iq c \ and J 4s , we only need to comment out the contraction 
terms (equivalently drop the terms proportional to B^) of Lf cl and Lf s of the "new" pro- 
grams. In this way, we found that the reconstructed IR subtraction terms from the "new" 
programs are identical with "old" ones for (a), (b), and (d). 
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However, the "old" IR subtraction term (c) I^i^ASm^ cannot be constructed by such a 
simple recipe from the "new" programs generated by GENCODEiV. Dropping the finite terms 



subtraction term l4b(i*)A5rri4 a from the scratch using the definitions of the fourth-order 
quantities I an d ASm^. Then, the result is compared with the integrand in the "old" 
program of AM°g d . 

C. I 4 ft( 1 *)A(5m4 a by the "old" 3-operation 

Let us first explain how J 4fe ( 1 *)A(5m4 a is obtained in the "old" program. In the old ap- 
proach, the IR subtraction term I^i^ASm^ originates from the 3i 2 37-operation. In addition 
to this term, J1237 operation yields the term / 4s AM 4(i . 

The IR-limit associated with the operator J1237 is given by 

z a + z b = 1 - 0(5), z 1 ,z 2 ,z 3 ,z 7 = 0(8), z 4 ,z 5 ,z 6 ,z c ,z d = 0(e), e ~ 5 2 , 5 — > 0. (53) 

In the neighborhood of this limit we have A 1 — 1 — 0(8), A 2 = 1 — 0(5), and V = 
0(5 2 ). As is discussed in the previous section, the result of 2?i237-operation includes both 
linear and logarithmic IR divergences. In particular, the linear divergence is associated with 
^137^1237 limit. If we apply the 3i373i237 operation, however, it subtracts the linear divergence 
correctly, but not the logarithmic divergence. Thus, we chose an ad hoc method in which 
the piece including linear divergence is separated out from the result of 3i 2 37-operation and 
put aside for a while. The remainder that contains only logarithmic divergence is named fk- 
The linear divergent piece is redefined so that it is defined on the subdiagram {1, 2, 3, 7, a, b} 
without decomposing it into two subdiagrams {2,6} and {1,3, 7, a}, which occurs in the 
naive ^137^1237 limit. This term is named The explicit forms of fk and /; in the old 



in M ^u^ASrriin is not enough to reproduce I^^ASm^a. Therefore, we reconstructed the 



FORTRAN program of AM 16 read: 




E + C + 5m 4a [f }+g s F 1 + Y 1 3(g s - V t )5m 4a [f ] + Y 





fl - ~n 




5m 4a [f ] 
U 2 V 4 



z 2 (l - A 2 f(-l + 6A 1 - 3Aj + 2A\) , 



(55) 
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where 



L 4s [F ] = (4Ai - 2A 2 )(1 -A ± + A 2 ,) + (-2 + A 1 A 2 ){1 - AA 1 + A\) 
F 1 = ^(2 - A 6 )/U + 2^6(1 - 2A 5 )/U + 5 56 (2 - A 4 )/C/ 

Fl = - Z 4 (B 45 (1 " 4>) + #46 + ,856^4)/^ 

+ z 5 (B 45 (l - A 6 ) - 4B 46 A 5 + B 56 (l - A 4 ))/U 

- z 6 (B 45 A 6 + B m + B m (l - A 4 ))/U 
V t = z 137 (l - Ai) - z 2 (l - A 2 ) 
E = 2A 4 A 5 A 6 - A 4 A 5 - A 4 A 6 - A 5 A 6 
Co = —3z c Zd/Us 

5m 4a [f ]=E + l-2A 5 

Y =z 4 (-A 4 + A 5 + A 6 + A 4 A 5 + A 4 A 6 - A 5 A 6 ) 

+ z 5 (l - A 4 A 5 + A 4 A 6 - A 5 A 6 + 2A 4 A 5 A 6 ) 

+ z 6 {A 4 + A 5 - A % - A 4 A 5 + A 4 A 6 + A 5 A 6 ) 
g s = z 4 A 4 + z 5 A 5 + zqAq . (56) 

The building blocks U,V, A, h By of the above integrands are obtained from those for M ie 
by taking the IR-limit associated with the 3i237-operation. Recall that in the IR-limit the 
subdiagram S consisting of the fermion lines 4, 5, 6 and photon lines c, d, and the reduced 
diagram Q jS consisting of the fermion lines 1, 2, 3, 7 and the photon lines a, b decouple from 
each other. Thus, the building blocks are actually the same as those obtained by taking the 
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UV-limit associated with the K456 limit. Their explicit forms are 



U = U S U g/S , 

Us = Z^QaiZc, + Z4 c ZQcii Ug/s = Zi 37a Z 2 b + Z 2 Z b 

V = V S + Vg/s, 

V S = ^456 - Z4A4 - Z 5 A 5 - ZqAq + X 2 Z cd , Vg/s = 21237 ~ ^137^1 ~ Z 2 A 2 + \ 2 Z ab , 

7 

A i = l-^z i B ij /U, i = l,---7, 

3=1 

B\\ = B ri = Byj = B 33 = B 37 = B 37 = z 2b Us, 

-B45 = Z 6 d Ug/s, B i6 = —Z 5 Ug/s, B 56 = Z ic Ug/s, 

Bij = for i e S and j e Q/S. (57) 

In the above z^... stands for z i± + z i2 + • • • and the electron mass is taken as a unit of 
mass scale (i.e., 1) and the photon mass is A. In the leading order of the J1237 limit L 4s [F ] 
tends to 4. The actual form of L 4s [F ] in Eq. (1561) was chosen so that the integral fk 
decouples into known lower-order quantities. This difference for /^[.Fo] is IR-finite. Note 
that 5m 4a [/ ] is related to the integrand of A5m 4a , namely, the UV- and IR-finite part of 
the mass renormalization constant <5m 4a given later in Eq. (|65|) . 

In order to clarify the structure of fk and //, let us split fk into the three parts fki, fk2, fm, 

f [ U \ 1 t rpi f Eg + Cp + gsF 1 + Yj t 3g s 8m 4a [f ] + Y 
fki = / (dz)g—L is {F ] — + 



AU 2 L J V V 3 V A 

fm = j (dz) g -^L 4s [F }5m, a [f Q } (± - ^I^zA^ j , 

fks = j (dz) g ^L^[F }8m, a [fo}^^^, (58) 
and fi into the two parts fa, fi 2 

f* = lf (dz)g S -^^z 2 A 2 (l - A 2 )(-l + 6A 1 - 3Aj + 2A\) , 

fa = -\J (dz)g 5 -^0^z 2 (l - A 2 ){2 - A 2 )(-l + 6A 1 - 3Aj + 2A\) . (59) 

The integrand fki was compared with the integrand I 456 M 16 generated by GENCODEiV by 
the spot-check method. We confirmed that the integral fki minus its K 2 limit, fk v , is equal 
to J 4s AM 4(l , which is listed as (b) of ( |52l) . It turns out that fk2 + fn is equal to I 4 6(i*) A£m 4a , 
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which is reconstructed from the lower order quantities in the next subsection. Thus, the 
difference between the "old" calculation and the reconstructed one is confined in /fe 3 + fa. 
It is IR-finite but contributes a nonzero value to AM,°' d . 



D. lMi*\A5m4 a reconstruction from the lower-order quantities 

In order to understand where this difference + fa came from, let us examine IR diver- 
gence structure of mass-inserted magnetic moment amplitude M 4 6(i*) in the "old" approach. 
i4&(i*) is defined from M 4 &(i*) as follows |f " 



(60) 



(61) 



i[36j,|38|: 

M 4 6 ( i*) = (1 - K 2 )(l - Ia«)(M46(i*) - M 4b (i*)[/]) + K 2 M Ab{1 , } 

+ l all {l - K 2 )M 46( i* ) [iVo + Z -f,E + Co] + M 4 6 ( i*)[/] 
= AM 4b(1 *) + (5m 2 M 2 ** + B^ V M 2 *) + / 46( i*) . 

The two terms in the second line of the r.h.s. define J 4 6(i*): 

hb(l*) = W 1 " K 2 )M 46(1 *)[iVo + Z - f, E + C ] + M 46(1 *)[/], 

where the function / is introduced in the first term in an ad hoc manner to subtract out 
the linear IR divergence coming from N + Z . The linear IR divergence is confined to the 
second term M 4 &(i*)[/], which has the form 

1 f f 

M 4 6(i*)[/] = -- / (dy)g /s 2 3 , 

J G/S Q/S 

f = -8y 2 A 2 {l - A 2 )(-l + 6Ax - ZA\ + 2A\) . (62) 
The explicit form of the first term of Eq. (1611) is 

l a uM 4b{1 4No + Zo-f,E + Co] = Jim f (dy) g/s ^^ [ -L- - Al) ) , (63) 

X ^°J u g/s \ lv g/s v g/s J 

but i where 

Vg/s = 2/137(1 - A,) + y 2 (l - A 2 ) + X\y a + y b ), 

(dy)g/s = dyidy 2 dy?,dy 7 dy a dy b 5{l - y 1237 - y ab ). (64) 
The finite part of the mass renormalization constant A5m 4a is 

A5m 4a = \ J (dy)s 5 ^, (65) 
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where <5m 4a [/o] is expressed in the same form in Eq. (156]) and 

(dy) s = dy 4 dy 5 dy 6 dy c dy d 5(l - 7/ 456 - y cd ). (66) 

Using the identity 

1 T(m + n) f\,f\, h , s t™- 1 ^- 1 

dt ds8{l-t-s)—— — ■ - (67) 



Ks% n /s r(m)r(n) j y v ^t^ + s^r' 

we can express the product of A5m 4a and the first term of i4&(i*) defined in Eq. f lBTl) in the 
same Feynman parameter space as that for the original amplitude M 16 

fH = A<5m 4a x I aH M 46(1 *)[iVo + Z - /, E + C Q ] 

L is [F ]Sm 4a [f ] ( 1 32i3 7 (l-Ai)\ 

4^ ^ V* J ' (68) 

which is identical with /fc 2 . Similarly, the contribution of the product of A<5m 4a M^(i*)[f] 
to M 16 is 

= I / Wo ^^ 2 (1 - A 2 )(-l + 6A! - 3Al + , (69) 

which is identical with fa. 

Therefore we find that the combination, —fk3 — fi2, is extra in the "old" AM°g d so that 
the correction term 

AMf ee 2(f k3 + f l2 ) 

= -2 x S y {dz) g 5 -^0t Z2 A 2 {l - Al f{l - A 2 ) , (70) 

where the overall factor 2 comes from time- reversal diagram, must be added to AM°g d . 
Evaluating it numerically, we obtain AMfg d = 0.029 437 8 (98), which is smaller than the 
current uncertainty of value B for M 16 in Table [TT] and cannot be detected by the direct 
comparison of value A and value B until the latter is evaluated more precisely. 

The difference between AMJ 1 ™ and AM°g d can be analyzed in the same manner. It is 
found that the difference is numerically not small for Mis' 

AM$ d = 2(l-K 5 )(f k3 + f l2 ) 

= -2 x 5 y {dz)g{ l - K 5 ) \^0fz 2 A 2 {l - A0 3 (1 - A 2 )} 

= -0.215 542 (19), (71) 
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where Ai, U, and V are denned in the J1237 limit of Mi§. Of course their explicit forms are 
different from those of M 16 . If we add AMf 8 dd to AMf 8 d , the value B for M 18 in Table [IT] be- 
comes 16.974 (21) and the difference between value A and value B is reduced to —0.006 (21), 
which is consistent with zero within the precision of numerical calculation. 

V. CONCLUSION 

The results described in this paper are summarized as follows: 

1) There was an inconsistency between the "old" integrals AM°g d and AM° 8 d and their 
residual renormalization terms. This inconsistency is resolved in this paper. 

2) Other 45 integrals of Group V of the "old" calculation are in good agreement with the 
"new" ones. 

3) Programs generated by GENCODEiV have no error for N = 8. Namely, the automation 
scheme has cleared the eighth-order test without difficulty 

The separation of the IR divergent and finite pieces in a given amplitude can be made 
arbitrarily. There is no overriding rule that dictates how to carry out such a separation. We 
only have to keep a record of what is subtracted as an IR subtraction term. All IR subtraction 
terms are summed up in the end and the arbitrariness in the choice of IR divergent part 
cancel out, leaving a finite contribution as a part of the residual renormalization. 

The important point is that the IR subtraction term prepared for the numerical calcu- 
lation and the one used to calculate the residual renormalization must be the same. What 
we found is that IWi*) used in the numerical calculation of and Mis and I^v) for the 
residual renormalization constant AM^mj had different forms in the FORTRAN programs 
of the "old" calculation. This is the reason why M 16 and M 18 had IR-finite but redundant 
contributions. 



The development of automatic code generator 1 15] was crucial in enabling us to 
discover the existence of extra IR subtraction terms in M 16 and M 18 on short notice. Adding 
the correction terms AMfg d and AMf dd to the "old" value, we find the entire contribution 
of Group V to be 




4 8) (GroupV) = -2.179 16 (343), 



(72) 
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which is in good agreement with the still preliminary value 

4 8) (GroupV) = -2.219 (53), (73) 

obtained by the "new" code generated by GENCODEiV. 

Due to the different forms of IR subtraction terms, the forms of the residual renormal- 
ization are also different in the "old" and "new" calculations. The residual renormalization 



terms in the "old" IR procedure are given by la, 136 . 
4 8) (GroupV) old = AM (8)old - 5AM (6)old A5 ; 
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- AA/f (4) {4AL (4) + 3A5 (4) - 9(AB 2 ) 2 } 

- M 2 {2AL (6) + A5 (6) - (10AL (4) + 6AB^)AB 2 + 5(A5 2 ) 3 } 

- AM^ASm^ 

- (M 2 * - M 2 4l}){A5m (6) - A5m w (5AB 2 + A5m 2 *)} 

+ M 2 A<W 4) (4AL 2 * + AB 2 * - B 2 * [I] ) , (74) 
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where 

47 



rold 



AM (8)oid = AM° 

i=01 
h 

AM (6)oid = AM 6 old 

x=a 

AM {4) = AM 4(i + AM 46 , 

AM< 4 *) = 2AM 4a(1 *) + AM 4a(2 *) + 2AM 46(1 *) + AM 46(2 *) 

ft 5 

AL {4) = 2AL 4c + AL Ax + 2AL 4s + AL 4i , 

h 

AB<® = J^^ASte, 

A5« = A£? 4a + A£ 46 , 
ft 

AcW 6) = ^r/ x A5m 6x , 

a;=a 

A^m*- 4 - 1 = A5m 4a + ASm^b 

1 for a; = a, 6, c, e, f , h 

(75) 

2 for x — d,g. 

The numerical values of the finite renormalization constants are listed in Tables IIVI and |V] 
and also in Appendix |A] B 2 * [I] is obtained from the 3 a ^-operation of the wave function 
renormalization constant B2*. 

The formula of the residual renormalization for the "new" calculation is much simpler 
than that for the "old" one. Since the mass renormalization is completed within the numer- 
ical calculation, the mass renormalization constant should not appear in the residual renor- 
malization. The exceptions are the vertex and wave-function renormalization constants that 
have self-energy subdiagrams. The mass inserted vertex (wave-function) renormalization 
constant L2*(I?2*) has no overall UV divergence. As a result, the i^-operation cannot pick 
up the renormalization terms proportional to L2*(I?2*)- It must be restored in the residual 
renormalization in order to carry out the complete on-shell renormalization. The residual 
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renormalization formula in the "new" approach is given by 



4 8) (GroupV) new = AM (8)new - 5AM (6)new AB 



2 



- AM^{3AL^ + - 9(AB 2 ) 3 } 

- M 2 {AL5 (6) - 6(AL (4) + AB^)AB 2 + 5(A5 2 ) 2 } 

+ M 2 A<W 4) (4AL 2 * + AB 2 , ) , (76) 



where 



ALB^ = AL (6) + Afi (6) + AL (4) AB 2 + A<W 4) 5 2 *[I], (77) 

AM (6)ncw = AM (6)old _ ^ _ M ^ [J]) A5m (4) _ M 2 AL (4) . (78) 

In Eq. (1761) the vertex renormalization constant AL^ and the wave function renormalization 
constant AB^ appear in the same weight for each order of the perturbation. It is because 
we have already subtracted one AL^ as an IR subtraction term. Calculating the combina- 
tion Ai>) + AB^ is much easier than calculating each of them separately. Because of the 
Ward-Takahashi-identity for the renormalization constants Lfo' + B^ n > = 0, many cancella- 
tions occur between two terms. Thus, we introduced a combined renormalization constant 
ALB^ 6 \ Its relation to the "old" renormalization constants AL^ and AB^ are given in 
Eq. (I77|) . More detailed definitions of ALB Sx for each diagram are given in Appendix IC 31 
The left-hand side of Eq. (1771) . ALB^\ was directly calculated with the p rog rams made by 



the automatic code generator for the residual renormalization constants 39j and obtained 

as 

h 

ALB^ = Vx&LB 6x = 0.100 86 (77) . (79) 

x=a 

This result was checked by comparing with the right-hand side of Eq. (ITT)) calculated using 
the residual renormalization constant for the "old" calculation. The sixth-order magnetic 
moment AM' 6 ' new was calculated with the programs generated by GENCODE./V and given as 

AM (6)now = Q.42610 (53). (80) 

See Appendix |A2] for the detail of AM( 6 ) new . 

As we have shown in the paper, the two results (1T2|) and (1T31) of the eighth-order contribu- 
tion from Group V diagrams are obtained by means of the totally independent calculations. 
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Further theory corrections to the eighth-order term of the electron g — 2 is very unlikely. 
The new theoretical prediction should be announced when we complete all the tenth-order 
calculation. 
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APPENDIX A: TEST OF GENCODEN BY LOWER-ORDER a e 

Although gencodeA was developed primarily to deal with the tenth-order g-type dia- 
grams, it can be readily applied to the calculation of fourth-, sixth-, and eighth-order g-type 
diagrams. Since these lower-order terms are also known from previous works, this serves for 
debugging of gencodeA. 

1. Fourth-order a e 

The fourth-order case is the simplest nontrivial example. The g — 2 receives correction 
at this order from four types of vertex diagrams, 4c, Ax, 4s and 4/ shown in Fig. H](b). 
Following the remark in Sec. Ill Al the sum of g—2 from the vertex diagrams 2M 4c ( s ) + M^rn 
is expressed as a quantity associated with a single self-energy diagram M 4a ( 4 f,) in Fig. 0] (a) 
via the Ward-Takahashi identity ( fl6i) . (The factor 2 is assigned to the diagram to account 
for the presence of the diagram which is related by reversing the orientation of the lepton 
line.) gencodeA creates two program sets for M 4a and within a few seconds on a 
generic Linux PC. These programs give the finite amplitudes AM 4(1 ( 4 h) as the sum of the 
unrenormalized Ward-Takahashi summed g — 2 amplitudes, also denoted by M 4a ( 4fe ) , and 
necessary UV and/or IR subtraction terms. The relation of AM 4a ( 4 b) to M 4a ( 4 b) is given in 
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M. 



4a 




M 



4b 




(a) 



(b) 



FIG. 4: Fourth-order q-type diagrams, (a) Self-energy- like diagrams M 4a and M 4b - (b) Vertex 
diagrams 4c, 4x,4s, and 4/. Their contributions to the magnetic moment are related to M 4a = 
2M 4c + M 4x , and M 4b = 2M 4s + M 4i . 



Appendix IB II It took about 10 minutes each to carry out their numerical integration by 
VEGAS 17| with ten million sampling points per iteration for 50 iterations on hp's Alpha 



machine. 

The values obtained in this way are 

AM 4a = 0.218 78 (35) , 
AM 4b = -0.187 73 (40) . 

The contribution from the fourth-order g-type diagrams is expressed as 

Af\q-type) = AM 4a + AM 4b - AB 2 M 2 



(Al) 



(A2) 



taking the residual renormalization into account. M 2 = 1/2 is the second-order correction 
to g—2 and AB 2 = 3/4 is the finite part of the second-order wave function renormalization 
constant B 2 . Substituting the numerical values (lAlj) for the formal expression f lA2j) . we 
obtain 



4 4) (g-type) = -0.343 95(53), 



which is in good agreement with the analytic value —0.344 166 • • 



(A3) 
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M 6c M 6d 







M 6e M 6f M 6g M 6h 



FIG. 5: Sixth-order g-type self-energy-like diagrams Mq x , x = a, . . . , h. The time reversal diagrams 
of M^d and M^ g are not shown here. The fermion lines of a diagram is named 1 to 5 from left to 
right. The vertex diagram obtained by inserting an external photon vertex into the fermion line i 
of the self-energy diagram 6x is named 6xi. 

2. Sixth-order a e 

The sixth-order diagrams can be evaluated in a similar manner and are found to give 



a result in good agreement with the numerically [20j and analytically [2l|] known values 
as follows. Fifty vertex diagrams of the sixth-order q-type diagrams are reduced to eight 
self-energy-like diagrams shown in Fig. [5] by means of the Ward-Takahashi identity and the 
time-reversal symmetry. It takes just one minute to create all eight FORTRAN programs 
for Mq x (x — a, b, . . . , h) by GENCODEiV on .hp's Alpha machine. Numerical evaluation was 
carried out on RIKEN's PC-cluster system (RSCC). After computation of 2 to 6 wall-clock 
hours with 16 Xeon-CPU's for each diagram to carry out a VEGAS integration with one 
hundred million sampling points per iteration for 450 iterations, we obtained 

h 

AM (6)ncw = ^ AM™ = 0.426 (11). (A4) 

x=a 

After continuing computation with one billion sampling points per iteration for 200 iterations 
for each diagram, we obtain the updated result AM' 6 '™ given in Eq. flHDl . The contribution 
of g-type diagrams to Af "* including the residual renormalization is given by 

4%-type) = AM (6)ncw - 3AM (4) A5 2 + M 2 {-A5 (4) - AL (4) + 2(A5 2 ) 2 } 

= 0.905 26 (53), (A5) 

where AB(L)^ is the sum of the finite parts of the fourth-order wave function (vertex) 
renormalization constants. See Appendix IC 21 for their definitions. The formula of the 
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residual renormalization (1A5I) can be obtained by using the definitions of the finite quantities 
AM 4(1 (j), AM 6x , etc., in Appendixes iBl and ICl The values of finite quantities such as AL 4c . . . 
are given in Table EJ The various finite pieces appearing in (1A5|) are 

A5 (4) = AB 4a + AB 4b = -0.437094(21) (A6) 
AL {4) = 2AL 4c + AL 4x + 2AL 4s + AL 4/ = 0.465 024 (17) (A7) 
AM (4) = AM 4a + AM Ab = 0.030 804 •• • (known exactly ) . ( A8) 

Eq. (1A5P again shows good agreement with the analytic result 0.904 979 • • • by Laporta and 
Remiddi [21 1. 



APPENDIX B: DIVERGENCE STRUCTURE OF THE MAGNETIC MOMENTS 

We briefly summarize our notation in Appendices [B] and O The relation between the 
unrenormalized amplitude M and the finite amplitude AM of magnetic moment is listed in 
this appendix. 

A symbol with a prefix A means a finite quantity. A renormalization constant with a 
superscript "UV", A uv , is the leading UV-divergent term of the on-shell renormalization 
constant A. A can be L, B, and 5m according to a vertex-, wave-function , and mass- 
renormalization constant, respectively. A VY is identical with A in Refs. jlq . [3^]. The 
subtraction terms proportional to a UV-renormalization term A vv are generated by the 
.fT-operations. A renormalization constant with a superscript "R" , such as L R sy is the resid- 
ual term defined in Eqs. ( 1271) and ( !29l . A of A R must be either L or 5m. B R is also 
defined accordingly, but there appears no B R term in the definition of a finite magnetic mo- 
ment amplitude. The subtraction terms involving A R are generated by the -R/i-subtraction 
operations. 

The subscript of M or A stands for the name of a diagram. A self-energy-like diagram 
of the second-, fourth-, and sixth-orders are called 2, 4a and 46, and 6a, 66, 6c, 6a", 6e, 6/, 6g, 
and 6h, respectively. (See Figs. H] and El) The 47 independent self-energy-like diagrams of 
the eighth-order are named from 01 to 47. (See Fig. [TJ) The fermion lines are always named 
from 1 to n — 1 from the left to right, where n is the order of the perturbation theory. 

The name of a vertex diagram is determined based on the Ward-Takahashi related self- 
energy diagram. When the vertex diagram is obtained by inserting the external photon 
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into the fermion line j of a self-energy diagram nx, this is called nxi. Thus, we have 
L>4xi, z = l,2,3, x = a, b, and L^ X i, i — 1, • • • 5, x = a, ■ ■ ■ h for the fourth-, and sixth-order 
vertex renormalization constants, respectively. The vertex renormalization constant of the 
second-order is named L 2 , since there is only one vertex diagam of the second order. In 
the early article 40], the fourth-order vertex diagrams are given other names. We follow 
the naming system Ref. 40j in this paper. The correspondence between two names of the 
fourth-order vertex diagrams is that 4al = 4a3 = 4c, 4a2 = 4x, 461 = 463 = 4s, and 
462 = 41. (See Fig. ij. 

When a two-point vertex is inserted into the fermion line i of a diagram nx, the resulting 
diagram is called nx(i*). Namely, * indicates the two-point vertex. 

The primed quantity, for example L 4x ^, is the derivative amplitude obtained by applying 
—Zi-J?- operation on the integrand, where Zi is the Feynman parameter assigned to the 
fermion line i. Note that L^') is equal to L 4x , but its UV-divergent part L^J^^ is not equal 
to 

^4x- Since second order quantities such as B2 and 5m,2 have only one electron line, it 
is not really necessary to distinguish different electron lines. We therefore use somewhat 
sloppy notations B2* and B21 instead of -62(1*) and B 2 iyy 

For L 2 , which contains electron lines 1 and 2, it is sometimes necessary to distinguish 
lines in which insertion is made. L 2 **t implies that two two-point vertices are inserted into 
the fermion line 1 of L 2 , while L 2 *t* means that one two-point vertex is inserted into the line 
1 and another into the line 2. M 4a contains three electron lines 1,2,3 and M 4a (i*+) means 
that two-point vertex insertion has been made twice in the electron line 1, and so on. 



1. Fourth-order magnetic moments 

The fourth-order magnetic moments are the same for both old and new approaches. The 
UV-finite amplitude are also given here. 

M 4a = AM 4a + 2 L 2 Y M 2 

AL,„ = A.U;„ 

M 4b = AM 4b + B^ y M 2 + 5m 2 M 2 + + L% M 2 
M 4b = AM 4b + Lf M 2 

In the new approach no explicit form of M 4 * is needed because the mass renormalization 
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is completed by the -R-subtraction operation. They are, however, listed here, since they are 
used in the old approach. 

M Aa (i*) = AM 4a (!*) + L 2 V M 2 * + h a (i*) 

M 4a (2*) = AM 4a ( 2 *) + ha(2*) 

M 4b( i*) = AM 46( i*) + (Sm 2 M 2 ** + B™ M 2 .) + I 4b{1 * } 

M Ab{2 *) = AM 46(2 *) + 8mi7 M 2 * + J 4b(2 *) + Lf M 2 * + M 2 *[I] 5m 2 * - 2 M 2 *[I] Lf 
2. Sixth-order magnetic moments by GENCODEiV 

The finite amplitudes of the sixth-order are given in the following. For simplicity, we 
drop the superscript "new" from AMg® w . 

M 6a = AM 6a + 2 5m 2 M 46(1 *) + 2 B™ M 4b - 5m 2 (5m 2 M 2 „ + B™ M 2 .) 

- B™ (5m 2 M 2 * + 5 2 uv M 2 ) 
+ 2 Lf s M 2 

M 6b = AM 6b + 5m 2 M 46(2 *) + B™ M 4b + Srrffi M 2 * + B™ M 2 - 5m 2 5m™ M 2 + 

- B™ (5m™ M 2 * + B™ M 2 ) 

+ Lf AM 4b + Lf t M 2 + M 2 ± A5m 4b 

M 6c = AM 6c + 2 L™ M 4b + 5m™ M 2 * + B™ M 2 - 2 L™ (5m 2 M 2 , + B™ M 2 ) 
+ Lf AM 4a + M 2 * A5m 4a 

M m = AM 6d + L™ M 2 + 5m 2 M 4a(1 * } + B™ M 4a + M 4b - B™ M 2 

- L 2 V (5m 2 M 2 * + 5 2 uv M 2 ) 
+ i£ M 2 

M 6e = AM 6e + 2 L u s v M 2 + <5m 2 M 4a(2 *) + B™ M 4a - 2 L 2 , v B™ M 2 
+ Ll M 2 

M 6/ = AM 6/ + 2 L™ M 2 + 2 L™ M 4a - 3 L 2 V M 2 

M 69 = AM 6s + L£ v M 2 + L™ M 2 + L™ M 4a - 2 L 2 V L 2 V M 2 

M 6h = AM 6/l + 2 M 2 



40 



3. eighth-order magnetic moments by GENCODEiV 



The finite amplitudes of the eighth-order are given in the following. For simplicity, we 
drop the superscript "new" from AM° ew . 

Moi = AM 01 + 2 Lf 1 M 6f + 2 L u c v M 4a + 2 M 2 - 3 (L" v ) 2 M 4a 

- 6 L 4 U C V M 2 + 4 (L^) 3 M 2 

M 02 = AM 02 + 5m 2 M 6/(1 *) + 5 2 uv M 6/ + L 2 V M M + L u s v M 4a + L u c v M 4fe + L 6 u v M 2 

- L 2 v (5m 2 M 4a( i*) + 5 2 uv M 4a ) - 5 2 uv M 4a - L 4c v (5m 2 M 2 * + 5 2 uv M 2 ) 

- fi 2 uv , M 2 - 2 L 4 U 7 M 2 - (L" v ) 2 M 4b + 2 L 2 V 5 2 uv M 2 
+ (L 2 V ) 2 (<5m 2 M 2 * + 5 2 uv M 2 ) 

+ ^/i M 2 

M 03 = AM 03 + 2 M 6d + 5m 2 M 6/(3 *) + 5 2 uv M 6/ + 2 M 2 

- 2 L 2 V (5m 2 M 4a(1 . } + B 2 UV M 4a ) - (L™) 2 M 4b - 2 L 2 V L u s v M 2 

- 2 £ 2 UV L^,,) M 2 + (L 2 V ) 2 (5m 2 M 2 * + £ 2 UV M 2 ) + 2 L 2 v 5 2 uv L 2 , v M 2 
+ L^ 3 M 2 

M 04 = AM 04 + 5m 2 M 6d( 3*) + £ 2 UV M M + 5m 2 M 6d(1 . } + 5 2 uv M 6d + L 2 V M 6a + M 2 

- 5m 2 (<5m 2 M 4a(1 ** } + 5 2 uv M 4a(1 * } ) - 5 2 uv (5m 2 M 4a(1 *) + 5 2 uv M 4a ) 

- 2 L 2 V (5m 2 M 46(1 * } + 5 2 uv M 4fe ) - 5 2 uv (L%, } + L^) M 2 

+ L 2 V 5m 2 (6m 2 M 2 ** + M 2 *) + L 2 V (5m 2 M 2 * + B^ v M 2 ) 
+ (5 2 uv ) 2 lg* M 2 
+ L%d3 M 2 + Lf dl M 2 

M 05 = AM 05 + M 6h + M 4a + Ma + L uv ^ _ 3 L uv L yv ^ 

M 06 = AM 06 + L 2 V M 6g + L 2 v M 6/ + L 4i v M 4a + M 2 + Ma 

- 2 (Dp 1 ) 2 M Aa - 2 M 2 - 3 L 4 u c v M 2 + 3 (L^) 3 M 2 
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M 07 = AM 07 + L™ M 6g + L™ M 6/ + L™ M 4a + L 6 U ^ M 2 + L™ M 2 
- 2 (L™f M 4a - 4 lUV M2 _ L uv L yv ^ + 3 (L uv )3 ^ 

M 08 = AM 08 + M 6c + 2 M 6d + 5m™ M 4a(1 * } + M 4a + L™ M 2 

2 (L™f M 4fe - (<5m u a v M 2 * + 5 4 u a v M 2 ) - 2 ( 5m2 + B uv M ^ 

- -uv r uv m _ B vv L uv Ma + 2 (L uv )2 (5m2 ^ + s uv Ms) 



- 2 L 2 v L 4s v M 2 - £> 4a ^ 2 , 
+ 2L^ £ 2 UV Lo3 v M 2 
+ AM 4a + Mml*) A5m 4a 



M 09 = AM 09 + <5m 2 M 6/(2 * } + 5 2 uv M 6/ + L 2 UV M 6e + L 4 U S V M 4a + M 2 + L 6 7 3 M 2 

-uv 

J 4c{2>) 



L™ (5m 2 M^**) + 5 2 uv M 4a ) - E 2 UV L^ v M 4a - 5 2 uv M 2 



B. 



UV rUV 



-^4c(l') ^2 



3 L™ L u s v M 2 + 3L 



™ 5 2 uv L 2 , v M 2 



+ -^6/ 2 ^2 



M10 = AM 10 + 5m 2 M 6d(2 *) + 5 2 uv M 6d + 5m™ M 4a(1 * } + 5 4 U 6 V M 4a 

+ L 2 V M 66 + L™ M 2 - 5m 2 5m™ M 4a(1 . } - B™ (5m^ v M 4a(1 * } + £# v M 4a ) 

- (5m 2 M 46(2 *) + £ 2 UV M 4b ) - £? 2 UV L u s | 2/) M 2 

- L™ (5m™ M 2 * + B™ M 2 ) - B™ L™ M 2 + L™ 5m 2 5m™ M 2 * 
+ L™ B™ (5m™ M 2 * + B™ M 2 ) + B™ B™ I™ M 2 

+ L* (AM Ab + L% M 2 ) - Lf M 2 + Lf d2 M 2 + M 4a(1 * } A5m 46 



M n = AM U + 2 5m 2 M 6d{5 * } + 2 B™ M 6d + 2 L™ M 4b 

- 5m 2 (5m 2 M 4a(1 * 3 *) + B™ M 4a(1 * } ) - B™ (5m 2 M 4a(1 * } + 5 2 uv M 4a ) 



2 L u s v (5m 2 M 2 * + B™ M 2 ) - 2 B™ L™ M Ab 
2 B™ L™ (5m 2 M 2 + + B™ M 2 ) 



'2 ^2 
+ 2 £ 6 d5 



M 2 
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Mi 2 = AM 12 + 2 5m 2 M 6a(n + 2 £ 2 UV M 6a + 5m 2 M 6a(3 *) + £ 2 UV M 6a 

- 2 5m 2 (<Sm 2 M46(i**) + £ 2 UV M 46(1 *)) - 2 5 2 uv (5m 2 M 46(1 *) + 5 2 uv M 4b ) 

- 5m 2 (5m 2 M 4fe(1 * 3 *) + M 46(1 *)) - B™ (5m 2 M 4b{1 * } + B™ M 4b ) 

+ (5m 2 ) 2 (5m 2 M 2 +++ + B™ M 2 ++) + 2 5m 2 B™ (5m 2 M 2 ++ + B™ M 2 +) 
+ (5 2 uv ) 2 (5m 2 M 2 * + 5 2 uv M 2 ) 
+ 2 Lf al M 2 + Lf a3 M 2 



Mi 3 = AM 13 + 5m 2 M 6fc( i*) + 5 2 uv M 6h + L™ M 4b + L™ M 2 



rUV 



(5m 2 M 2 * + £ 2 UV M 2 ) - 5 2 uv M 2 



M 14 = AM 14 + 5m 2 M 6g(5 *) + E 2 UV M 69 + L 2 UV M 6(i + L 4 U Z V M 4b + L 6 u rf v 3 M 2 



uv 



•uv 



■uv 



L 2 UV (5m 2 M 4a(1 . } + 5 2 uv M 4a ) - L 4 U , V (5m 2 M 2 * + 5 2 uv M 2 ) - £? 2 UV M 2 
(L^) 2 M 4b - L 2 V L 4 U S V M 2 + (L^) 2 (5m 2 M 2 * + 5 2 uv M 2 ) + L 2 V 5 2 uv I™ M 2 



+ ^5 ^2 



M 15 = AM 15 + 5m 2 M 6fl(1 *) + S 2 UV M 6s + L 2 UV M M + L 4 U C V M 4fe + L 6 u a v 2 M 2 



uv 



-uv 



•uv 



- L 2 V (<5m 2 M 4a(1 * } + B 2 UV M 4a ) - L 4 U C V (<5m 2 M 2 * + £ 2 UV M 2 ) - 5 2 uv M 2 



>uv 



•uv 



uv 



>UV rUV 



(L^) 2 M 4b - L. 



UV rUV 



Lt: m 2 

uv 



+ (L 2 UV ) 2 (5m 2 M 2 , + 5 2 uv M 2 ) + L 2 UV B 
+ ^1 M 2 



-UV dUV rUV 



L 2 U , V M 2 



M 16 = AM 16 + 5m 2 M 6c(1 * } + £ 2 UV M 6c + <5m u a v M 46(1 * } + £ 4 u a v M 4b + 2 L 2 V M, 

- 5m 2 (5m™ M 2 „ + B™ M 2 *) - B™ (5m™ M 2 , + B™ M 2 ) 

- 4 L™ (<5m 2 M 46(1 . } + E 2 UV M 4fc ) 

+ 2 L 2 V 5m 2 (5m 2 M 2 ++ + B™ M 2 .) + 2 5 2 uv L 2 V (5m 2 M 2 * + 5™ M 2 ) 



6 a 



+ M 2 + Lf s AM 4a + M 



46(1 * } A<5m 4a - A<5m 4a M 2 
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Mir = AM 17 + <5m 2 M 6e(1 . } + 5 2 uv M 6e + <5m 2 M 6(i(4 *) + £ 2 UV M 6d + L^J M 4b 
+ M 2 - 5m 2 (<5m 2 M 4a(1 * 2 * } + 5 2 uv M 4a(1 *)) 

- £ 2 UV (5m 2 M 4a(2 * } + £ 2 UV M 4a ) - L u s v (<5m 2 M 2 * + E 2 UV M 2 ) - 2 5 2 uv I™ v) M 2 

- 5 2 uv M Ab + 5m 2 5 2 uv M 2 * + (5 2 uv ) 2 (L™ + L"X) M 2 
+ M 2 + L* 4 M 2 

Mis = AM 18 + 5m 2 M 66(1 *) + 5 2 uv M 6b + <5m 2 M 6a(2 *) + 5 2 uv M 6a + <5m u 6 v M 46(1 *) 

+ B 4 U 6 V M 4b - <5m 2 (<5m 2 M 4b(1 * 2 *) + 5 2 uv M 46(1 * } ) - E 2 UV (<5m 2 M 46(2 *) + £ 2 UV M 4fe ) 

- 5m 2 (5m u 6 v M 2 ** + fl™ M 2 *) - 5 2 uv (<5m u 6 v M 2 * + B™ M 2 ) 

- 5m 2 <Jm£Y M 46(1 *) - 5 2 uv (tfm^ M 46(1 *) + M 46 ) 

+ <5m 2 8m%y (5m 2 M 2 ** + B^ Y M 2 *) + B^ Y 5m^ (5m 2 M 2 ** + B^ Y M 2 *) 
+ B 2 UV (5m 2 M 2 * + 5 2 uv M 2 ) 

+ L* 61 M 2 + Lf s (AM 4fe + M 2 ) + L* 2 M 2 - M 2 
+ (M46(i*) - L f* M 2) A5m 4{ , 



M19 = AM 19 + 2 M 2 



M 20 = AM 20 + L 2 Y M 6h + I%J 2 M 2 + M 2 - 2 L 2 V M 2 



M 21 = AM 21 + 2 M 2 



M 22 = AM 22 + M 6g + L u c v M 4a + L™ M 2 + L^ 2 M 2 

- (LD 2 M 4a - 3 L 2 V L 4 U C V M 2 - m 2 + 2 (L^) 3 M 2 

M 23 = AM 23 + 5m 2 M 6ft(2 *) + 5 2 uv M 6/l + L u e2 M 2 

UV ji^ D UV / rUV , rUV 



+ L 6d4 M 2 - B 2 (L 4x(ll) + L 4x(2 ,)) M 2 

+ ^2 M 2 



M 24 = AM 24 + 5m 2 M 6g{2 * } + B 2 UV M 6 , + M 4a + M 2 
+ L™ M 2 - B^ M 4a - 5 2 uv M 2 
— 2 £ 4s v L 2 Y M 2 — B 2 Y ^ 4( ^3/) Af 2 + 2 L 2 v B 2 Y L 2 Y M 2 
+ L%2 M 2 
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M 25 = AM 25 + 2 Lf M 6g + 2 L$J 2 M 2 



(L^)' M 4a - 2 L 4l 



UV rUV M 



2 L 



UV rUV 



L 4 U C V M 2 + 2 (L^) 3 M 2 



M 26 = AM 26 + 2 M 6c + 2 L u c v M Ab + 5m^ v M 2 * + E 6 U / M 2 - 3 (L^) 2 M 4b 

- 2 (5m u a v M 2 * + 5 4 u a v M 2 ) - 2 L u c v (5m 2 M 2 * + 5 2 uv M 2 ) 
+ 3 (L^) 2 (<5m 2 M 2 * + £ 2 UV M 2 ) 

+ AM 6/ + M 2 * A5m 6/ 

M 27 = AM 27 + 5m 2 M 6gm + 5 2 uv M 6 , + L uv M 6e + L u d v M 2 + L^ 2 M 2 

- 5 2 uv L™, , M 2 - 5 2 uv M 2 - 2 L u s v M 2 - L uv (5m 2 M 4a(2 . } 



+ £? 2 UV M 4a ) + 2 L 2 UV £ 2 UV L 2 U , V M : 
+ ^4 M 2 



M 28 = AM 28 + Sm 2 M 6c(2 * } + 5 2 uv M 6c + L 2 UV M 6b + L 4 U / M 4b + 6m£ M 2 * + M 2 



- L 2 UV (5m 2 M 46(2 *) + 5 2 uv M 4b ) - 5 2 uv L 2 U , V M ib - 5m 2 Sm^ M 2 * 

, .u 2 , - ^ (5m 4fe v M 2 * + M 2 ) 



^ UV M 2 * + S 4 U J 1M M 2 ) - (5m 4 u fe v M 2 * + S?v 



J V"" t 4a(l') J "2* -r ^ 4a (i') ^"2J ~ -^ 2 

(<5m 2 M 2 * + M 2 ) + Lf 1 5m 2 5m^ M 2 + 



L uv 



+ L 2 v 5 2 uv (5m 2 u , v M 2 , + £ 2 U , V M 2 ) + £ 2 UV L 2 U , V (<5m 2 M 2 * + 5 2 uv M 2 ) 



,uv 



>uv 



UV rUV 



? uv 



+ L* 2 M 2 - L 2 Lf c M 2 + L 2 (AM 6d + Lf c M 2 ) + M 2 * A5m 6d 



M, 



= AM 29 + 2 5m 2 M 6e(2 *) + 2 5 2 uv M 6e + 2 L u a \ M 2 - 5m 2 (5m 2 M 4a(2 **) + 5 2 uv M 4a(2 *)) 
B 2 UV (5m 2 M 4a(2 .) + £ 2 UV M 4a ) - 2 5 2 uv (L^, + Lgf iw) ) M 2 + 2 (5 2 uv ) 2 M 2 



+ 2 L* 2 M 2 



M 30 = AM 30 + 2 <5m 2 M 6b(2 *) + 2 £ 2 UV M 6fe + <5m u a v M 2 * + E 6 u a v M 2 

- 5m 2 (5m 2 M 46(2 **) + 5 2 uv M 4b(2 *)) - £ 2 UV (<5m 2 M 46(2 *) + 5 2 uv M 46 ) 

- 2 5m 2 (JmSfi*) M 2 * - 2 5 2 uv ((Jm^j M 2 * + BY h J v) M 2 ) 
+ 2 5m 2 5 2 uv <Sm™ M 2 * + (5 2 uv ) 2 (Sm^ M 2 * + B^,7 M 2 ) 

+ 2 L^ fe2 M 2 - 2 ^ 4s M 2 + (AM 6a + 2 M 2 ) + M 2 * A5m 6a 
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M 3 i = AM 31 + 2 M 2 



M 32 = AM 32 + M 2 + L" v 2 M 2 



M 33 = AM33 + 2 Lgv M2 



M 34 = AM34 + L V J M Aa + I™ M 2 + M 2 - 2 M 2 



M 35 = AM35 + M 6h + M 2 + Ma _ 2 L uv L uv ^ 



M 36 = AM 36 + L uv M 6 , + L 4 7 M 4a + M 2 + L 6 u g v 5 M 2 - (L 



uv 



uv 



uv 



UV\2 



4a 



- 3 M 2 - L 2 UV L 4 U C V M 2 + 2 (L 2 uv f M 2 



UV rUV 



UV\3 



M 37 = AM 37 + 2 M2 



M 38 = AM 38 + 2 M 4b + 5m^ M 2 * + M 2 - 2 L^ v (<5m 2 M 2 * + 5 2 uv M 2 ) 



,uv 



uv 



-uv 



>uv 



+ Lf AM 6h + M 2 * A5m 



M 39 = AM 39 + L 2 UV M 6g + L 4 U C V M 4a + Lj£ M 2 + L 6 U ^ M 2 - (L 2 v ) M 4a 
- 3 L 4 U C V M 2 - L 2 V L 4 U , V M 2 + 2 (L^) 3 M 2 



M 40 = AM 40 + L 2 V M 6c + M 4fe + L 4c v M 46 + 5m" v M 2 * + 5 6 u g v M 2 

- 2 (L^) 2 M 46 - L 2 V (<fm™ M 2 * + 5 4 u a v M 2 ) - (5m 2 M 2 * + 5 2 uv M 2 ) 

- L 4c v (5m 2 M 2 * + 5 2 uv M 2 ) + 2 {L^f (5m 2 M 2 * + 5 2 uv M 2 ) 
+ AM 69 + M 2 * A5m 6g 



M 4 i = AM 41 + 2 L 2 V M 6e + <Jm^ M 4a(2 * } + 5 4 u a v M 4a + 2 L 6 u v M 2 



- 2 L 2 UV (5m 2 M 4a(2 * } + 5 2 uv M 4a ) - 4 L 2 UV L 4 U / M 2 - 2 5 4 u a v L 2 U , V M 2 
+ L 4x AM 4a + M 4a(2 .) A5m 4a 
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M 42 = AM 42 + 2 L™ Me6 + 5m uv M4b(2i) + B vw M ^ + Jm uv M ^ + ^uv M2 

- 2 (5m 2 M 4b(2 *) + B™ M 46 ) - 2 (5m™ M 2 , + £ 4 U 6 V M 2 ) 

- 5m™ 5m™ M 2 * - B™ (5m™ M 2 * + B™ M 2 ) 

+ 2 L™ <5m 2 5m™ M 2 * + 2 L^ v B™ (5m™ M 2 * + B™ M 2 ) 

+ Ll AM 4a - (Lf) 2 AM 4a + Lf (AM 6c + Lf AM 4a + M 2 * A5m 4a ) 

+ M 2 , (A5m 6c + Lf 5mf a ) + M 4fr(y) A5m 4a 

- t^ra^ A5m 4a M 2 * - A5m 4a Lf M 2 * 



M 43 = AM 43 + 5m 2 M 6h(3 .) + £ 2 UV M 6h + 2 L 6 u e v 2 M 2 - 2 E 2 UV L™, } M 2 



M 44 = AM 44 + 5m 2 M 69(3 *) + 5 2 uv M 6fl + L 4 U S V M 4a + M 2 + L^\ M 2 

-UV A/f R UV rUV 
v 42(2') iW 2 — -t>2 ^4c(2') 



fi? v L^; V M 4(7 , — i?y v L™oi\ M 2 — B™ L™ (0 ,\ M 2 — 2 L u s v L™ M 2 



'2 ^2' J "4a ±-"2 

+ 2 5 2 uv L™ M, 
+ L%s M 2 



M 45 = AM 45 + 5m 2 M 6c(3 *) + B™ M 6c + 2 L™ M 4b + 5m™ M 2 * + B™ M 2 
- 2 B™ L™ M 4b - 5m 2 5m™ 2 ^ M 2 * - B™ (5m™ 2l) M 2 , + B™ 2f) M 2 ) 



- 2 L u s v (5m 2 M 2 * + B™ M 2 ) + 2 B™ L™ (5m 2 M 2 * + B 2 
+ L* 3 M 2 + Lf (AM fe + M 2 ) + M 2 * A<5m 6e - L* Lf M 2 



M 2 ) 



Vf 46 + 5m 2 M 6e(3 *) + 5 2 uv M 6e + 5m™ M 4a(2 *) + B™ M 4a 



2 L™ M 2 



M 46 = AM, 

- 5m 2 5m™ M 4a(2 *) - B™ (5m™ M 4a(2 *) + B™ M 4a ) - 2 B™ L™ 2I) M 2 

- 2 B™ L™ M 2 + 2 B™ B™ L™ M 2 

+ L* 3 M 2 - M 2 + (AM 4fe + Lf M 2 ) + M 4a(2 * } A5m 4b 
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M 47 = AM 47 + 5m 2 M 6K3 *) + v3 2 uv M 66 + 8m^ M 4b(2 *) + 5 4 U 6 V M 46 + 5m^ M 2 * 
+ B™ M 2 - 5m 2 5m™ M 4fe(2 *) - B™ (5m™ M 4fe(2 *) + B™ M 4b ) 

- 5m 2 5m™ 2 ^ M 2 * - B™ (5m™ {2/) M 2 , + B™ 2 , } M 2 ) - 8m™ 5m™ M 2 * 

- B™ (5m™ M 2 * + B™ M 2 ) + 5m 2 (5m™) 2 M 2 * + B™ 5m™ 5m™ M 2 * 
+ B™ B™ (5m™ M 2 * + B™ M 2 ) 

+ M 2 * [A5m 6b + Lf {5mf b - (5m 2 8m* + B™ 5m%,)}\ 

+ (AM 6b + M 2 * A5m 4b + AM 46 + L% M 2 ) + M 4 R b(2i) A5m Ab 

+ (AM 4b + L% M 2 ) + Lf b3 M 2 - A5m 4b 8m 2 * M 2 * - L% A5m 4b M 2 * 

- (Lf) 2 (AM 4b + Lf M 2 ) - 2 Ll L% M 2 + (L*) 3 M 2 

APPENDIX C: DIVERGENCE STRUCTURE OF THE RENORMALIZATION 
CONSTANTS 

1. Second-order renormalization constants 

L 2 = L™ + L 2 , = L 2 = I 2 

B 2 = B™ + B 2 , B^ = B 2 = —I 2 + A B 2 

4~ B^ — AB 2 
B 2 * = —2 L 2 * , L 2 * = I 2 * + AL 2 * 
B 2 ** = — 2 (2 Zv 2 **t + L 2 *t*) 
5m 2 * = 5m™ + I 2 + A5m 2 * 

2. Fourth-order renormalization constants 

Lix = L™ + J 4x + AL 4x 

L ic = L™ + J 4c + AL 4c + L™ L 2 

B^a = B™ — I 4x + AB 4a + 2 L 2 V B 2 — 2 7 4c 
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Lai = L\f + hi + (Lf) 2 + AL 4l + L 2 L 



2 



Las = L V J + I 4s + AL 4s + 5m 2 L 2 * + £ 2 UV L 2 , 

B 4b = B 4b Y + AB41, + 5m 2 B 2 * + B 2 V B 2 > + Lf B 2 — 2 / 4s — J 4 ; 



ALB 4a = 2 + Lf x + Bf a = 2 AL 4c + AL 4x + AB 4a 

ALB 4b = 2 L* + + Bf a - Lf AB 2 = 2 AL 4s + AL 4l + AB 4a 



ALB& = ALB 4a + ALB 4b = AL (4) + A£< 4 > 



3. Sixth-order renormalization constants 



L 6al = Lg al + 2 5m 2 L 4s( i*) + 2 L 4s(r) - 5m 2 (5m 2 L 2 **t + £ 2 '*) 

- 5 2 uv {5m 2 L 2 „ + 5 2 uv Z 2 „) + 

Lgo.2 = Lf a2 + L 2 V L 4s + <5m 2 I^i*) + L 4 ;(i') — L 2 Y (5m 2 L 2 * + L 2 <) + 

-^6a3 = L% a3 + 2 (5m 2 L 4s( i*) + B 2 V L 4s{V) ) - <5m 2 (5m 2 L 2 * t * + L 2 '*) 

- £? 2 UV (5m 2 L 2 ,* + 5 2 uv L 2 „) + 

^661 = + ^2 -^4s(2*) + B 2 V L 4s ( 2 /) + 5m u b v L 2 * + 5 46 v L 2 > - 5m 2 5m 2 * v L 2 * 

- t> 2 [dm 2 , L 2 * + ±S 2 , L 2 i) + L m 

L§b2 = Lf b2 + 5m 2 L 4 /( 2 *) + 5^ L 4 /( 2 /) + L 4s v L 2 — B 2 V L 2 ^ L 2 + L^ b \ 

-^663 — -^663 + Lj 2 -^4/ + ^Al Lj 2 ~ \^ 2 ) ^2 + ^ 

L 6cl = Lf cl + 2 Z 4s + 5m u a v L 2 * + B V J L 2 , - 2 L 2 v (Sm 2 L>. + ^ v L,) + A, 

-^6c2 — -^6c2 + Lj 2 -^4/ + -^ 4c ^2 ~ \^ 2 ) ^2 + Ll , 

T - fR i fUV 7 , rUV 

^6c3 — -^6c3 + -^4a; ^2 + -^6c3 

^6di = L^ dl + <5m 2 L 4c (x*) + B 2 V L 4c (i/) + L 2 V L 4s — L 2 V (Sm 2 L 2 * + B 2 V L 2 >) + Lg d \ 



-uv 

J 6b3 



UV 
6cl 



■UV 
7 6c2 
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T — T R _i_ 7~ uv 7 _|_ rUV 7 frUV\2 7 J- t 
^6d2 — J^6d2 + L 2 -^4c + -^2 -^4/ _ \J^2 ) ^2 + U 



UV 
6d2 



L 6d3 = Lf d3 + 5m 2 L 4c (i*) + £ 2 UV L 4c( i + L 4s - (<5m 2 L 2 * + 5 2 uv L 2 ,) + L 



uv 

6d3 



L&dA — Lfd4 + ^ m 2 ^4a;(l*) + -B^ ^4x(l') + -^6d4 



-^6cZ5 — -^6d5 + ^ m 2 -^4c(3*) + -B^ ^4c(3') + £4/ ^2 ~ -B^ -^2' V ^2 + ^6d5 



-^6el — -^6el + <5m 2 -^4c(2*) + -B^ -^4c(2') + I/ 2 — -B^ L 2 + L 



UV 
6el 



-^6e2 — Lf e 2 + ^ m 2 -^4x(2*) + -B^ ^ix{2') + -^6e2 



r _ rR I rUV r , rUV 

-^6e3 — -^6e3 + -^2 -^4x + -^6e3 



UV 



r _rRirUVr , rUV r / rUV\2 r i r 
-^6/1 — -^6/1 + -^2 -^4c + -^4c -^2 — ^2 J ^2 + ^6/1 



r _ rR i rUV 7 , rUV 
-^6/2 — -^6/2 + -^2 ^^x + -^6/2 



UV 



r _rR lOrUVr /rUV\2 r , r 

-^6/3 — -^6/3 + 1 L 2 L *c — {^2 ) L 2 + -^6/3 



^6gl — -^6ol + -^2 ^4c + -^4 C ^2 ~ \J^2 ) -^2 + 



UV 
6gl 



t _ rR _i_ r UV 

-^6g2 — -^ 6 g2 + L %g2 



t _ rR _i_ t UV 

-t>6 9 3 - ^ 6 <?3 + ^693 



r _rRirUV7 , rUV 
-^6g4 — -^6g4 + -^2 + -^6 S 4 



r _rR 1 rUV 7 irUVT frUV\2 7 , r 
-^6g5 — -^6g5 + -^2 -^4c + ^Al ^2 ~ \^2 ) 2 ~r ^ 



UV 
695 



r _rR , rUV 7 1 r UV 



r — rR 1 rUV 

^&h2 — ^%h2 + ^6h2 
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uv 

6h3 



B 6a = B* + 2 (5m 2 S46(i*) + 5 2 uv B 4b(v) ) - 5m 2 (5m 2 B 2 ** + 5 2 uv £ 2 '*) 

- 5 2 uv (5m 2 5 2 « + 5 2 uv £ 2 /,) + £ 6 u a v 

£? 66 = Bf b + 5m 2 £ 46(2 *) + 5 2 uv B^a/) + 5m u , v E 2 * + B™ B 2 , - 5m 2 8m%? B 2 * 

- 5 2 uv (5m™ B 2 * + B™ B 2I ) + B™ 

B 6c = Bf c + 2 B ib + 5m™ B 2 * + B™ B v - 2 L™ (5m 2 B 2 * + B™ B 2 ,) + 5 6 u c v 

= Bf d + 5m 2 -B 4a (i*) + B 2 y i?4a(i') + Bib + L™ B 2 — L 2 V (<5m 2 _B 2 * + S 2 B 2 >) 

R uv r uv n , R uv 

- t> 2 L 2 , t> 2 + U 6d 

Be e = B^ e + 5m 2 -B 4a ( 2 *) + -B^ £ 4a ( 2 ') + 2 £ u s v B 2 — 2 B 2 V L 2 ? B 2 + _B^ e v 

R _ RR i O rUV r _,_ 9 rUV p o (tUV\2 g , pUV 
±>6/ — -Dgf + Z -L 2 -D 4a + ^ -^4c ^ 2 — 6 [L 2 J ti 2 + iJgj 

R _ dR , rUV S , rUV p , rUV 5 9 /rUV\2 p , rUV 
-D6g — -D6 g + ^2 -°4a + -^4c ^2 + -^4/ ^ 2 - I [L 2 j & 2 + £> 6g 

R dR 1 o rUV r 1 rUV 



The residual renormalization constants ALB Qx for each diagram appearing in Eq. (ITT)) 
are defined in the following equations. 



AL£? 6a 




.R 

'6a 



- 2 L 4 A£? 2 



ALB 6b - 2 L 6bl + 2 L 6fe2 + L 6b3 + B ( 
- Lf (Bf b + 2 Lf s + L^) - L 



,R 

66 



A5 2 + (L^) 2 A5 2 
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ALB 6c = 2 Lg cl + 2 L^ c2 + Lg c3 + Bf c 

- Lf (Bf a + 2LI + Ll) 

ALB ed = Lg dl + Lf d2 + Lg d3 + L^ M + L^ d5 + 5^ 

- Ll AB 2 

ALB 6e = 2 L^ el + 2 Lg e2 + Lf e3 + Bf e 

- Ll AB 2 

ALB 6 f = 2 L^ fl + 2 Lg /2 + Lg /3 + Bl f 
ALB 6g = Lg 9l + Lg 32 + Lg g3 + Lg ff4 + Lg g5 + 



ALB 6h — 2 Lg M + 2 Lg\ 2 + Lg\ 3 + Bf h 

h 

ALB {&) =J2vxALB 6x = AL (6) + A5 (6) + AL (4) A5 2 + AcW 4) £ 2 *[7] 
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